|hep-th/ 9803162 
IPM-98-284 
March 1998 



Distributed Systems of Intersecting 
Branes at Arbitrary Angles 

R. Abbaspur* t[] and H. Arfaei*-tQ 

* Institute for Studies in Theoretical Physics and Mathematics, 
P.O. Box 19395-5531, Tehran, Iran. 
^Department of Physics, Sharif University of Technology, 
P. 0. Box 19365-9161, Tehran, Iran. 



Abstract 

A 'reduced' action formulation for a general class of the supergravity solutions, 
corresponding to the 'marginally' bound 'distributed' systems of various types of 
branes at arbitrary angles, is developed. It turns out that all the information re- 
garding the classical features of such solutions is encoded in a first order Lagrangian 
(the 'reduced' Lagrangian) corresponding to the desired geometry of branes. The 
marginal solution for a system of iV such distributions (for various distribution 
functions) span an iV dimensional submanifold of the fields' configuration (target) 
space, parametrised by a set of N independent harmonic functions on the transverse 
space. This submanifold, which we call it as the 'i/-surface', is a null surface with 
respect to a metric on the configuration space, which is defined by the reduced La- 
grangian. The equations of motion then transform to a set of equations describing 
the embedding of a null geodesic surface in this space, which is identified as the H- 
surface. Using these facts, we present a very simple derivation of the conventional 
orthogonal solutions together with their intersection rules. Then a new solution for 
a (distributed) pair of p-branes at SU(2) angles in D dimensions is derived. 
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Introduction 

BPS configurations of intersecting branes have been the basis of many recent develope- 
ments in string and M-theory. This is mainly because of their essential role in estab- 
lishing string dualities, counting the entropy of extremal black holes and constructing 
super symmetric gauge theories (for reviews see 0-0 and referrences therein). The con- 
struction of supergravity solutions corresponding to such configurations of branes have 
been implemented using four main approaches. These include: supersymmetry tech- 
niques j|, H H |7|, [15], [l6j , duality transformations || ^ [K| [15], [16[] , dimensional reduction 



and oxidation |TT], [T3|, [15], and lastly direct solving of the bosonic field equations of 



supergravity fTT|, [18], [19|, ^T], |22], If we are going to study, in a unified manner, a 
wide class of these solutions in various types of D-dimensional supergravities, containing 
all types of p-branes , the only systematic approach is to use the direct method. In this 
approach one considers a purely bosonic theory consisting of two sectors, the 'gravita- 
tional sector' including a metric tensor and several dilatons, and the 'form-fields sector' 
consisting of antisymmetric tensors ( 'form-fields' ) of various degrees. Such a model can 
be imagined to be the bosonic sector of a supergravity theory in D dimensions, obtained 
from the 1 ID-theory via the processes of dimensional reduction and truncation [|TTJ. In 
the usual approach, one writes the complete set of the field equations for the corresponding 
.D-dimensional action, and tries to solve them by inserting a set of ansatze appropriate 
to the given configuration of p-branes . In spite of its logical completeness and relia- 
bility, such an approach involves extra complications due to the appearance of coupled 
Einstein-'form-fields' equations. Typically, introducing an ansatz with a set of field equa- 
tions, means 'reducing' the 'actual' degrees of freedom (including spacetime dependences) 
to the 'relevant' (i.e. excited) ones describing the physical situation of interest. So an 
alternative approach is to form a 'reduced action' describing only dynamics of the excited 
degrees of freedom, and try to solve their equations of motion. However, in this manner we 
will lose the equations of motion corresponding to 'irrelevant' degrees of freedom, which 
appear as the constraints between relevant ones. For this reason, not every solution of the 
'reduced theory' is also a solution of the original one, but it is easy to see that the converse 
is true. The situation is different when certain 'subspaces' of the 'configuration space' are 
concerned^. It turns out that supergravity BPS saturated solutions describe two types of 
brane systems || 0: z)marginal configurations and ii)non-marginal configurations. In the 
first class the binding energy of the bound system of branes for their arbitrary separations 
vanishes so that they don't feel any total force. In the second class the binding energy 

3 This is comparable to a similar fact in the differential geometry that: Every geodesic line of a 
Riemannian space, lieing on an arbitrary surface, is also a geodesic line of the surface itself, while the 
converse is not necessarily true. However for certain subspaces called 'geodesic surfaces' the inverse 
theorem is also true. 
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is a negative quantity and the total force between each pairs of branes is attractive. As 
a result the bound state can be stable only at zero separations and the corresponding 
solutions are spherically symmetric around the common centre of branes in their trans- 
verse space. The marginality condition in the first class, imposes heavy constraints on the 
form of the corresponding solutions. They must satisfy certain 'extremality' and 'no-force' 
conditions. Further, all the field variables must be written totally as functions of a set of 
'independent' harmonic functions | |2~4fl . It will be seen in this paper that, the (first order) 
'reduced Lagrangian' of the theory £, on the 'subspace' H of such solutions, identically 
vanishes. This property, with some insights from the differential geometry, enforces the 
idea that every solution of the reduced theory on the 'if -surface' must describe a solution 
of the original theory as well. It seems that all the physical information regarding these 
solutions comes essentially through the corresponding reduced actions. In this paper we 
try to elaborate this idea, to the extent that it can be used for determining the solutions 
corresponding to the orthogonal, as well as the non-orthogonal configurations of p-branes . 
We will study solutions with uniform distributions of branes along their relative transverse 
directions, but with arbitrary distributions along their overall transverse directions. We 
shall refer to such systems as the 'distributed' systems, which are described in terms of a 
set of density functions in this paper. Continuous distributions in a low-energy model, can 
be interpreted as the long-distance limits of periodic (or un-periodic) arrays of branes with 
relative separations of order \fa! in the high-energy string theory. In the low-energy limit 
the Kaluza-Klein modes corresponding to the compactified relative transfers coordinates, 
are averaged and the solutions will depend only on the overall transverse coordinates |TJ . 



We begin our study of these solutions in section 1, by a reformulation of the problem of 
two different orthogonal branes, in the framework of a reduced theory. We will introduce 
a set of constraints and insist on their determining role in solving the field equations of 
the reduced theory. We will find in this approach, in agreement with those of [H]], that 
the field equations are completely replaced with a set of algebraic constraints determining 
parameters of the solution. The relations between couplings and dimensions emerge as 
the consistency conditions of these algebraic constraints. As a byproduct, a Diophantine 
equation governing the allowed marginal intersections of super p-branes is obtained. It 
is argued that the usual relation between the dilaton couplings and the corresponding 



form-fields degrees in supergravity [|l7j] , is a direct result of supersymmetry. In section 2 
we generalize the ideas of section 1 to a system of multiply intersecting orthogonal branes. 
It will be seen that the algebraic constraints in matrix representation have a universal 
form. Then the solution for these constraints together with their consistency conditions in 
compact forms will be given. In section 3 we lay the foundations of a more general theory 
of the marginal brane solutions including systems of branes at angles. General ansatze for 
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the form-fields describing purely electric-type branes, in terms of the ansatz for the metric 
are derived. It is shown that these ansatze are dependent on a set of 'structure constants' 
, which encode the geometry (i.e. angles) of the brane system. A general formula for these 
constants, in terms of the asymptotic form of the metric is derived. The concepts of the 
'if-surface' and an 'if-basis' are defined. A general formula for the reduced Lagrangian 
(RL), describing such generalized systems is derived. Using this RL the equations for 
the embedding of the 'if-surface' in the configuration space of the field variables will be 
written. Heuristic derivations for the generalized versions of the constraints in sections 1 
& 2 are given and the role of the C = equation is emphasized. In section 4 we apply 
the general framework of section 3, to the problem of arbitrary branes at angles. Then 
explicit solutions for the marginal configuration: p H p = (p — 2) using this formulation is 
derived. During the way of this derivation we encounter an integrability condition, from 
which the relation between angles (besides other information) is obtained. In section 5 a 
formulation of the no-force conditions appropriate to the framework of section 3 will be 
given. It is shown that these conditions in combination with the extremality conditions 
gives rise to a class of constraints, derived earlier on an ad hoc basis. The various (mass- 
less) fields contributions to the static long-range potential (between a brane system and a 
brane probe) are separated. The consistency conditions of section 1 emerge as the balance 
conditions among the long-range forces. In section 6 the formulas for masses and charges 
are presented and shown, as is expected that, these become proportional to the integrals 
of the density functions. In Appendix A a general first order Lagrangian for gravity is 
formulated and its application for the derivation of the RL's is indicated. In Appendix 
B we present two methods for classifying the solutions of the Diophantine equation for 
intersections for arbitrary spacetime dimension D, and then give the classifications for 
D = 4,6,10,11. In Appendix C the angles between a pair of branes in terms of the 
asymptotic form of the spacetime metric are defined. Finally in Appendix D we give the 
proofs of some if -surface identities. In the following sections, we will use a model theory 
with an action of the form: 

I D = j d D x (y=gR - \iy V f - \ ± e^J^j (0.1) 

where jF r 's denote {p r + 2)-form field-strengths defined as: T r = dA r , and = 
(p +2)\ FrM 1 ...M Pr+2 Fr 4l ''' Mpr+2 ■ Upon fixing the non- vanishing form-fields and the respective 
couplings a r , we are able to handle various supergravity models with suitable truncations. 
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1 Two Orthogonal Branes 



To show the essence of the reduced theory formulation, we begin our discussion by con- 
centrating on a simple example. We consider a simple distributed system constructed on 
the basis of a pair of orthogonal [] {d\ — 1, d 2 — l)-branes intersecting (overlapping) f] over 
{5 — 1) dimensions in the spacetime of arbitrary dimension D. The (minimum) dimension 
of a homogeneous hyperplane filled with the uniform (parallel) distributions of the two 
branes obviously is 

d = d 1 + d 2 -5 (1.1) 

The overall transverse space has D — d = d + 2 dimensions. We will use a set of Cartesian 
coordinates (x 1 *, y™, y™\ z a ) as defined in table (1). 



coordinate 


dimension 


tangent (v) 




5 


v || d\ , v | d 2 


V? 


8\ = di — S 


v || d\ , v J_ 5 


V2 


5 2 = d 2 -5 


v || d 2 , v J_ 8 


z a 


d + 2 


v _L di , v _L d 2 



table{l) 



By construction all the components of any physical field for such a configuration will 
depend only on the coordinates (z a ) of the transverse space. We can also allow a 'trans- 
verse' distribution of branes, which unlike the 'longitudinal' distributions may be inhomo- 
geneous. Throughout this paper we will describe such 'transverse' distributions using the 
density functions p\(z), where A labels^ various distributions in the problem (e.g. A = 1, 2 
for now). It must be remarked here that such distributions are possible only when each 
pair of the parallel constituent branes can form stable bound states with zero binding 
energy. It is widely believed that this 'marginal' binding property, is a characteristic 
of supersymmetric bound states P, |5|, [T|, p7|. However it will become clear during this 
section that the distributions of non-supersymmetric p-branes in the non-supersymmetric 
theories are also possible. Also, we will see that bound states of various such distributions 
become possible whenever they have a suitable number of common directions. 
To construct field theoretic classical solutions, as usual a set of ansatze for the metric and 
the form fields are required. The ansatz for the metric is fixed by means of its isome- 
tries defined by the set of commuting Killing vectors: (d^,d m ,d m '), and its isotropies: 
SO(5 — 1, 1) x SO(5i) x SO(5 2 ) x SO(d + 2). Here we have demanded Lorentz invariance 

4 Clearly the parallel situation d\ C d-i is a special case where 8 = d\. 

5 In this paper we use the term 'intersecting branes' in equal footing with 'overlapping branes', i.e. 
different branes may have separated centres in their transverse space. We use the world— volume notation: 
d\ H di = S for indicating such intersections. 

6 In this paper we will use the indices A, A', etc. and k, k', etc. exclusively for the distribution quantities. 
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along the common directions. This requirement is believed to be a characteristic of 



super symmetric solutions [17, H. However the solutions discussed here are not limited 



to supersymmetric theories. Hence we write the ansatz for the metric as 

ds 2 = e 2Bo ^dx^d Xfl + e 2Bl{z) dy™dy™ + e 2B ^ z) dyf dyf + e 2C{z) dz a dz a (1.2) 

where dx^dx^ = r] flu dx fJj dx 1 ' and 77^ = diag(—l, +1, +1). We should fix the ansatze for 
the form fields as well. For this purpose we must specify that whether the (d\ — 1, d 2 — 1)- 
branes carry electric or magnetic charges, i.e. whether they couple to [d\ + 1, d 2 + l)-form 
field strengths or to their dual (d± + 1, d 2 + l)-forms. To be specific, we assume here both 
of them to be of the 'electric' type (generalizations to the cases of magnetic- or mixed-type 
branes are then straightforward JT7J). Then, as long as d\ 7^ d 2l the Lorentz & rotational 
invariance along the subspaces fix the (d\, d 2 )-form potentials to 

Ai = ±e Xl ^{dx^)A(dyT) 

A 2 = ±e X2{z \dx> l )A(dyf) (1.3) 

where the ± signs discriminate between branes and anti-branes, and the differentials 
in the parentheses are abbreviations for the wedge products over the indicated indices. 
Whenever d\ = d 2 and simultaneously oti = a 2 , the above invariances fix the ansatze 
for Ai and A 2 up to any linear combinations of the above two expressions. Because in 
this case we have no mean for distinguishing the forms in our model. Therefore, we can 
replace these two potentials by a single one written as: A = ±Ai ± A 2 . (Note that this 
property is restricted to the orthogonal branes, where we have: T 2 = T\ + !F 2 , so that 
sum of the kinetic terms of (A\,A 2 ) are replaceable by that of A only.) 

However when d\ = d 2 but a\ ^ a 2 , such a replacement is not possible. For example 
in Type IIB theory, the bound state of orthogonal D- and NS-strings excites the 2-form 
fields of the RR and NSNS sectors: A 2 and B 2 respectively. Despite being orthogonal, 
these 2-forms can not be traded for a single 2-form, as we have: a a = -«b = 1/2 FH , TD| . 
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The reduced Lagrangian (RL) formulation 

To form a reduced action describing the physical situation of interest, first of all we have 
to specify the relevant degrees of freedom. In the case of interest these are X\ with 
A = 1, 2, and 4> a = (B , B\, B 2 , C, cp) with a — 1, 5. Here we have unified the dilaton 
with the independent metric variables into a column vector <p a for simplifying our later 
formulations. The RL for (<ft a (z), X x (z)) is then found by inserting the ansatze (|1.2|) 
and ( |1.3|) in the action (p.l|) , and integrating over the relative transverse coordinates: 
(x^, y™, y™')- This, using the formula (|A.11| ), gives the RL 



C = C 



G 



where Cq and £p, respectively the gravitational (including the dilaton) 
and form field parts of the RL, have the expressions 



Cf 



l/2e 2G tt a pt 
1/2 E 



D 2(X X ~F X ) 



A=l,2 



(ox x y 

s defined by 



;i.5) 



Here d = d a and (G, Fa) are linear functions of 

G{4>) = g a r 

Fx(ct>) = fx«r 

and (Q a f3,9a, f\a) are matrices with constant entries depending on the dimensions, 



1.6) 



aj3 



( 6(5-1) 

55! 

55 2 
d(d+ 1) 
\ ' 



86i 

6 l (d + l) 




Qa = (5/2 SJ2 5 2 /2 d/2 
/i a = (5/2 5 x /2 -5 2 /2 -d/2 0) 
/2« = (5/2 Si/2 d 2 /2 -d/2 o; 

The equations of motion resulting from this RL are 
5C 



55 2 d(d + 1) 

5x5 2 8 x {d + l) 

5 2 (5 2 -l) 5 2 (d + l) 

5 2 (d+l) d(d+l) 



0) 



o \ 





-1/2/ 



;i.7) 



5X X 
5C 



d(e 



°dX, 



D 2(X A -F A ) 



d(e 2G n aP d^) - g a e 2U n 



2G ( 



/3 7 ( 



(dX x ) 



+ E / 

A=l,2 



2(X X -F K ) 



(dX x 



(U 



Obviously a general solution for such a set of nonlinear 2 nd order equations is not easy to 
find. However we will see that a class of solutions, corresponding to a pair of localized 
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distributions p x (z), can be found using simple tricks. The suitable boundary conditions 
for these configurations are those of the asymptotic flat metric, constant dilaton and 
vanishing field strengths at infinity, i.e. 

4> a \z^oo = 4> a o ) Xx\z-kxi — X\ (1.9) 

In what follows we always take Q O = 0, which is always possible by choosing suitable 
scales of coordinates (for a = 1,...,4), and of the form fields (for a = 5) (Note that a 
constant shift in the dilaton field can beabsorbed by scaling of the form-fields, as far as 
the Chern-Simons terms are not introduced in the action ( p.l|) ). Then of course we will 
not be at liberty on the definitions of X Ao 's, since as will be seen below, they have absolute 
physical meanings in terms of the mass to charge ratios. In particular when X Ao = we 
will find that our solution describes bound state of super p-branes. 



Solving the equations using the constraints 

Although the equations (|1.8|) seem formidable to solve, if the exponentials can be taken 



(consistently) to be constants, then they will be greatly simplified. So we look for the 
solutions satisfying the constraints 

G{(j)) = 

X A -F A (0)=X Ao (1.10) 



The constants on the right sides have been fixed using the boundary conditions (|1.9|) with 
<ft a = 0. Later we will find that these constraints have the physical interpretations as the 
'extremality' and the 'no-force' conditions respectively. Using only the latter constraints 
in the X A 's equation of motion we find that 

d 2 e - Xx = (1.11) 

(More precisely equation ( |1 . 1 1| ) must be written as a Poisson equation: d 2 e~ Xx ^ = 
—K\p\{z). However, assuming that distributions fill only a local region of the z space, eq. 
( |1.11| ) can still be used in the outer 'empty' region.) 

The general solution to this equation with the boundary condition as in ( |1.9| ) is written 

as 

C -**M = K ?H x (z) (1.12) 
where the constants k x and the harmonic functions H x are defined as 

K X = e Xxo (1.13) 

H x {z) = 1 + J d^ 2 z' Px {z')G^ 2 {z, z') (1.14) 



7 



Gj +2 (z, z') is the Green's function of the (d + 2)-dimensional (flat space) Laplacian, 



G d+2 (z,z') = l \*-*'\ ^ ' (1.15) 

I — m | z — z | , d = 

The masses and charges of these distributions are found in section 6 to be proportional 
to the integrals of corresponding p x (z)'s, showing that these describe 'densities' of the 
distributions. Applying our constraints to the second equation in (|1.8| ), we find 



Qap&tP - gjl^d^.df + £ K X 2 f Xa (dX x ) 2 = (1.16) 

A=l,2 

We can solve these equations for a 's using a further ansatz which is proposed naturally 
by our constraints ( |1.10| ). Linearity of these constraints in both (<p a ,X x ) suggests that 
we take the linear ansatz: 

4> a {X) = - £ Cx(X x -X X0 ), (1.17) 

A=l,2 

with suitable choices of the coefficients £™, we will find a simultaneous solution of the 
constraints ( TT7T0 ) and the equations of motion ( |1.16| ) for <fi a (z). Assuming X x s to be 
independent (which is always possible by independent choices of the distributions p x (z)), 
our constraints (|1.10|) and equations of motion (|1.16| ) are respectively replaced by the two 
sets of algebraic equations: 

g-£ A = o 

fx-£x> = -6xx> (1-18) 

and 

Ca^Ca'^0 , A ^ A' (1.19) 

where we have used the matrix notation for simplicity. These equations have a solution 
for (£« g) as, 

j, a 2/ ^1 _ ^1 ^1 ^1 \ 

Kx ~ Kl [ ~2D , ~2D , 2D , 2D ,ai) 

fa 2r d 2 d 2 d 2 d 2 . , . 

= K2( -2D'2D'-2D'2D' a2) (L20) 

where D = D — 2 and = D — d x . In addition they give a set of 'consistency 
conditions' , 

2 1 d\d\ 
ai = 2D 
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2 1 G?2<^2 

ai« 2 = ^-*) (1-21) 

which express in a different notation f\ the result of [^TJ for the equations governing the 
intersections (usually called the 'intersection rules' ) of the two electric type branes. They 
are derived here simply using a set of constraints. Generalisations of the relations ( |1.21|) 
can be obtained in a variety of models with various scalars and form-fields correspond- 
ing to the various types of intersecting branes E7I]. We will see in section 5 that 



equations (|1.21|) have simple interpretations in terms of no-force conditions among the 
same/different types of branes. For now let us concentrate on the first two of them, i.e. the 
equation 

„ 2 1 d X d X n 00 , 

ax = ^ ~ ~2D (L22) 
This relation is not specific to the bound states of orthogonal distributions; even the 
solution to a single distribution requires such a condition. The simple interpretation 



of ( |l-22j) is obtained by calculating the mass M x and charge Q x of a (d x — l)-branes 



distribution (formulas (|6.11| ) and (|6.13|) ), 

M A = n x 2 j d* +2 z Px {z) 

Qx = n x jd d+2 zp x (z) (1.23) 
We observe that the physical meaning of K\, 

* - £ (1-24) 

VA 

is the mass to charge ratio. This justifies our previous statement that X Xo = In k x can not 
be set by hand to arbitrary values. Using this ratio in (|1.22|) , we find a relation between 
the three 'charges' f\oi a (d\ — l)-brane: 

Qx 2 = a x 2 M x 2 + ^M A 2 (1.25) 

This relation will become physically meaningful if we recall that the long range forces 
between two parallel identical (d x — l)-branes have the forms 

M x a x 2 M x 2 Q x 2 



7 In and most of the other references, the cha's are twice the a x s in this paper. Also k x s are often 
traded for another set of parameters defined by: A x = 4/k\ 2 . 

8 These 'charges' include the gravitational, dilatonic and form-field charges, which are proportional to 
M\, a x M\ and Q\ respectively. Compare to jl(| for a similar discussion in the context of black hole 
solutions. 
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whereFc, Fd and Fp represent the gravitational, dilatonic and (d\ — l)-form forces and r 
is the transverse distance between two branes. In this way the equation Ql.25 ) is nothing 



but the simple statement that the gravitational and dilatonic attractions of two (d\ — li- 
braries must be cancelled against their (d\ — l)-form repulsion, to allow stable bound 
state. Alternatively we might consider a single (d\ — l)-brane instead of a pair of them 
(i.e. Px(z) ~ 5 dx+1 (z)). In that case again we would recover equation (|1.25|) , but of 



course with a different interpretation as the balance between the various long-range forces 
between different parts of a single brane. In other words, if a brane has enough 'charge' , 
its form-field repulsion will prevent its collapsing due to the gravitational/dilatonic 'self 
attractions' (when it has been slightly 'bended' ). 

For the third relation in (|1.21|) , a similar argument must yield an interpretation in 
terms of the pair-wise no-force conditions of the intersecting branes (see section 5 and the 
discussion of fT0|l). 

An interesting special case of the above solution concerns the bound states of super 
p-branes. In such a case the unbroken 1/2 supersymmetry of each individual (d\ — 1)- 
brane causes it to have the minimum mass:MA = \Q\\ (i.e. k,\ = ±1). Such branes 
are of special interest, because in addition to being classically stable against the long- 
range forces (which is guaranteed by ( |1.25[ )), they are also stable against loop corrections 
and Hawking radiations |39], ffO]| . Strictly speaking such (supersymmetric) objects occur 



(only) in supersymmetric theories like the supergravity and superstring theories. So not 
surprisingly supersymmetry of the theory requires special relations among its parameters 
so as to render it compatible with the existence of such extended objects. For the model 
at hand a relation of this kind, restricting the couplings and dimensions, is obtained by 



putting K\ = 1 in equation ( 1.22[) . Thus in agreement with fl7|l , we find that 



«a 2 = 1 - d -§ (1.27) 

This equation reproduces (up to a sign) the dilaton couplings to the form-fields in all 
D — 10 supergravity theories with one dilaton field. In D = 11 supergravity which has 
no dilaton field (a a = 0), hence we find that only the dual 2 and 5-branes (corresponding 
to 3 and 6-form potentials) are possible as supersymmetric extended objects. 
Returning to our double brane system, from ( 1.21|) we obtain, 



(2D - d 1 d 1 )(2D - d 2 d 2 ) = (5D - d x d 2 ) 2 (1.28) 

We shall refer to this Diophantine equation as the 'intersection rule' in this paper, since it 
specifies the suitable dimension of the intersection subspace for marginal configurations. 
In Appendix B methods for classifying the solutions of ( |1.28| ) for arbitrary D are proposed 
and the possible solutions for D = 4,6, 10, 11 are indicated. It must be noted however that 



10 



there exist also BPS non- marginal configurations [|| 0, which do not obey the intersection 
rule ( |1.28 ). As has been indicated in ||, the non- marginal solutions in D = 10 can be 



constructed from the marginal ones using T- and S-duality transformations. For theories 
in D < 10 and D = 11 this can be done using the methods of dimensional reduction 
and oxidation fLlf . However a general theory of such solutions for arbitrary spacetime 
dimensions is presently absent. 

Summary of the solutions for orthogonal di R d<i = 5 system 

The solutions for the metric and dilaton are found from ( |1.17|) , referring to the definitions 



of the variables and values of the parameters in ( |1.20| ). Also the form-potentials are 
already known using their expressions by ( |1.3|) and (|1.12|) . The results for the general 
case (arbitrary k\) are 

ds 2 = Hf dllT) Hf d2lT) {H^ 2 R^ 2 dx^dx^ + H{ Kl2 dy™dy™ + H^dyfdyf + dz a dz a } 
A 1 = KrH^idx") A (dy?) 

A 2 = K 2 H 2 \dx^) A {dyf) (1.29) 

Therefore we see that the solutions of the reduced theory, satisfying suitable constraints, 
are in complete agreement with those of plj p3fl obtained by solving 'directly' the field 



equations of the original action ( |0.1| ). As pointed in the introduction, this feature is a 
general property of the marginal solutions satisfying the C = condition, which is seen 
to be manifest for the above solution. 



2 Generalization to the multiply intersecting branes 

To illustrate the generality and simplicity provided by applying the methods introduced 
in section 1, we will apply them to obtain the solution for iV (orthogonal) intersecting 
branes. The steps will be exactly parallel to those of section 1. We will use, for simplicity, 
of a model with one dilaton field <p and iV form-potentials A\ of degrees d\ with A = 
1,...,N, which are coupled to iV orthogonal branes of electric type with the respective 
world-volumes dimensions. Again the whole of the branes are assumed to have uniform 
distributions within a (^-dimensional world-volume and arbitrary distributions transverse 
to it. Clearly the (minimum possible) dimension of this wold- volume (d) is 

d = E rf A - J2 + E rf AA'A» - - (2-1) 

A A<A' A<A'<A" 
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where A, A', ... = 1, N, and d\\>, g?aa'A"'s, etc. are the dimensions of double, triple, etc. 
intersections respectively. We may assume that the time direction is included in all the 
d\s so that the configuration is static, though it is not an essential restriction. We will use 
the Cartesian coordinates (x l ) and (z a ) for the homogeneous subspace and the subspace 
transverse to it respectively. 



The ansatze for ds 2 and Ax's 
We choose a set of world- volume coordinates (x % ), whose directions span the various 
constituent branes' world- volumes. In these coordinates we can write 



ds 2 = e 2Bi{z) Vij dx i dx j + e 2C(z) dz a dz a 
A x = ±e Xx(z) /\. J dx l 

where r/ij = diag(—l, +1, +1) ; i,j = 0, d — 1. 



(2.2) 



The RL and all that... 

The relevant degrees of freedom are 4> a = (B i: C : (p) and X\. The corresponding RL is 
written as 

C = \l1e 2G ^^d^.d^ + 1/2 ]T e 2{x *- F ^\dX x ) 2 (2.3) 



A=l 



where as in section 1 we have defined (N + 1) linear functions of (0 a ) as 

G(0) = g a r 



and the (d + 2)-dimensional constant matrices as 



/ o 

1 



a/3 



1 





1 1 

d+1 d+1 
V 



1 
1 





d+1 
d+1 

d+1 



d+1 d(d+l) 




\ 




-1/2/ 



9a 

The symbol e\i is 



1/2(1, ..,1, d, 0) 

1/2(±1, ±1, -d, -2a x ) = l/2(e Ai , -d, -2a x ) 



+ 1 
-1 



i E d\ 
i £ d x 



(2.4) 



(2.5) 



(2.6) 
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It is important to note that as far as we have not included any extra term (e.g. Chern- 
Simons terms) in our action (p.l|) , the general form of the corresponding RL for orthog- 
onally intersecting branes always resembles that of (2.3). The only difference of such 
models appears in the definitions of the variables 4> a and of the constants (f Xa ,g a , ^ap)- 
Hence a marginal orthogonal solution in every case is obtained, by following the proce- 
dure introduced in the previous section. In particular one obtains the same equations of 
motion, constraints, harmonicity condition for e _XA 's, ansatz for (f> a (X), and consistency 
equations as the equations flTgP , flTToD , ([Till) , flTlT]) and ( jOg) & flEISp respectively 



(but now with A = 1 , . . . , N) . The final link in this chain always is to solve the set of 
algebraic equations ( |1.1S| ) & ( |1.19| ) simultaneously for finding (if exists) the consistent 
values of our £°'s. Of course the existence of such a consistent solution, often will be 
accompanied by a set of relations between the parameters of the model itself, which as we 
have seen, lead to physical information about the underlying theory. The solution to the 
equations (|1.18|) & (|1.19|) together with the necessary consistency conditions are found by 



simple algebraic manipulations to be, 

£ x = K X 2 n- 1 (s + f x ) (2.7) 



gn-^g + fA) = o 



K\K\> 



{21 



In particular for the model at hand, replacing (fl, g, fx) with those of (|2.5| ) one obtains 
from O 

(2.9) 



Zx = ^ 2 + 2D' 2D' ax) 



where we have defined 



The first consistency condition in (p 



a x a X ' - 
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i e d x 
, i £ d x 

becomes identity 
1 / d x d y 



(2.10) 



while the second gives 



K X K X > 



D 



XX' 



(2.11) 



This is a closed form of the three equations in (|1.21|) for each pair of (d x — l,d\/ — l) branes 
intersecting over (d X y — 1) dimensions. This equation has the simple interpretation that: 
for iV branes to be marginally bounded, the various 'two-body' forces among them must 



vanish separately. If for example d XX i X » had been appeared in (|2.11|) , we should have 
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interpreted it as the vanishing of the 'three-body' forces which obviously was weaker than 
the above condition. The final expressions for the solutions are written as 



A' 



N 



ds 2 



n h x k > dx/D {(n h x -^ ^^dx 3 + dz a dz a } 

A=l A=l 



N 



A=l 



Ay 



'X12) 



where H x s are N independent harmonic functions. This reproduces the result of [21 



m a 



simple way, thereby gives another proof of the 'harmonic functions rule' |p4| . Again k x s 
measure the mass to charge ratios. So when = ±1 we obtain the solution for bound 
state of super p-branes . In this case equation ( |1.28| ) generalizes to a pair-wise intersection 
rule of the form 



{2D - d X d X ){2D - dydy) = {d xx ,D - d x d> 



(A ^ A') 



(2.13) 



a relation restricting different dimensions involved. 



3 General framework for (marginal) multiply inter- 
secting branes 

Before proceeding to discuss the case of branes at angles we need to somehow develop the 
tools of the previous sections to include the more general cases of marginal bound states of 
p-branes. The key point about the corresponding solutions is that they can be described 
in terms of a set of independent harmonic functions, in one-to-one correspondence to 
the transverse distributions of branes. (by definition each distribution consists of only 
the same type parallel branes). In order to give a general formulation, we consider a 
spacetime of arbitrary dimension D containing a homogeneous (but in general anisotropic 
) time-like surface of dimension d < D — 3. This subspace can be the world-volume 
of various types of p-branes with arbitrary relative orientations (and/or velocities) and 
separations, which are filled (and/or moving) uniformly within a (d — l)-dimensional 
space-like volume. Such a spacetime admits a set of Killing vectors which generate the 
symmetry group: R d x SO{d + 2), where the two factors represent the translational and 
rotational invariances along the distribution space and the transverse space respectively. 
We can decentralize arbitrarily the positions of different uniform distributions along the 
transverse space to obtain non- uniform distributions. This will break the SO(d + 2) 
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rotational symmetry, but not the R d translational invariance. So as in section 2, we can 
introduce two sets of Cartesian coordinates (x l ) and (z a ), along the distribution and the 
transverse space respectively. In these coordinates the ansatze for the metric and various 
form-potentials generally are written as 

ds 2 = h ij {z)dx i dx j + e 2C(z) dz a dz a 

A r = Ah(z)(da*) (3.1) 

where in the last formula f\ A r is a (p r + l)-form potential, generated by the set of all 
p r -branes of the same type, which are assumed to carry 'electric' charges of the corre- 
sponding form-fields. This choice makes all the ^4 r 's to have non-vanishing components 
only within the d dimensional world-volume ( Compare to a static charge distribution in 
electrodynamics). Generalization to systems involving branes with 'magnetic' charges is 
straightforward. 

To construct the solution for a system of N marginally bound distributions, we have to 
take all the field variables, say A 's, as functions of N independent harmonic functions: 
H X ;X = 1,...,N, 

<P A (z) = <^{H x {z\ H N (z)) = <P a (H(z)) (3.2) 

Using this in the equations of motion for (j) (;z)'s, we will obtain a new set of equations 
describing (j) A (H)'s as functions on the space of harmonic functions. We will use this 
technique partially throughout this paper to solve the equations of motion . A systematic 
approach to this viewpoint, and its further consequences will be presented in a forthcom- 
ing paper Q . 



Generation of the suitable ansatze for A r "'s 
First of all, we can use the above principle to obtain general ansatze for the form- 
potentials, which in all the non-orthogonal cases has a non-trivial expression. For this 
purpose we need the kinetic term of A r in the RL, which is calculated, using the ansatze 

(0 

C F = -±<i«»***{h i *)dAfo.dAfo (3.3) 

where G is defined through 

e 2G = 4^~ge- 2C = Vhe ic (3.4) 



9 In this formula we have introduced a convention, according to which the set of indices of a p-form 
field, as well as the wedge products over them, are indicated by a representative within a parenthesis. In 
this notation, a summation over all the repeated indices within a parenthesis, including a division by p\ 
is implied. 
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In this formula and later on, we use the convention introduced in (3.1). So in (|3.3|) a multiplication 



over /i l:P s followed by summations over (i) and (J), and a normalization factor: l/(p r + 1)! are implied. 
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and h = —det(hij). So we find the equation of motion for A r as 

d( e 2G+2a ^(h ij )dA r U) ) = (3.5) 
which using dAfy = dH\d x A r ^ and 

d 2 H x = (3.6) 

becomes 

d xl (e 2G+2a ^(h ij )d x A r {j) )dH x .dH y = (3.7) 

The left hand side of ( |3.7|) is written as a quadratic polynomial function of the variables 
dH x . Assuming independence of H x s, this equation implies 

d x {e 2G+2a ^d x ,A r {j) ) + d x ,{e 2G+2a ^d x A r {j) ) = (3.8) 

which upon integration yields, 

F r x {r) = e 2G+2 ^(h^d x A r U) = + c r $H x> (3.9) 

where cjf* and c$ are two sets of constants with the latter having the antisymmetry 
property 

c$ + c$ = (3.10) 

Now we argue that each of the A r, s components must be dependent only on those subset 
of H x s describing the charge distributions of the corresponding branes , i.e. those with 
d X = Pr + 1- In other words, the form-field associated to the charges of a given type 
must be independent of the charges of the other types. This is a consequence of the fact 
that the equation of motion for A r , i.e. d * (e 2oLr ' p J rr ) = *J r contains only the currents 
of the (corresponding) p r -branes as the source term. This feature, according to ( |3.9|) . 
corresponds to the existence of a consistent solution with the properties 

for d x ^p r + l J 9 ^ = r ° _ (3.11) 

[ C A — C AA' — U 

For d\ = p r + 1 we see from ( fj.lU| ) and (|3.11[) that c^y in ( |3.9| ) only mixes H x s of the 
corresponding A r sector. In fact consistency requires that all c^'s in each sector to be 
vanishing. So we can write (|3.9|), after inverting it for d x A r ^, as 

d x Afo = e-* a ^<$ i) {h ij ) (3.12) 

This gives a general electric ansatz for T T in terms of the dilaton and the metric compo- 
nents. This is simply achieved by taking the curl of A r in (|3.1|) using (|3.6|) and ( |3.12|) , 
the result is 

F = ufaAidx*) (3.13) 
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where ojLSs are a set of 1-forms defined by 



= dA r {i) 



-2(G+a r (f>) r(i) 



(3.14) 



If we know the i^A-dependences of ((p,hij,C), then we can integrate ( |3.12j ) to express 
-4[i)'s (and so ^4 r 's) as explicit functions of (H\). Of course this is possible only whenever 
certain integrability conditions are satisfied. These conditions are written simply as the 
Bianchi identities: dT r = which are equivalent to 



duj, 







(3.15) 



Equation ( |3.15| ) leads to a set of first order partial differential equations for the unknown 
functions (<p(H), C(H), hij(H)). Since these functions must also satisfy the corresponding 
equations of motion , the above construction (in principle) determines the possible rela- 
tions among s. This in turn specifies the allowed set of angles between the various 
p-branes, since as we will see shortly, all these constants have expressions in terms of 
angles. 



Some special cases 

Before proceeding further, it worths to examine here the above result for some well known 
marginal BPS solutions. For simplicity we will consider solutions with only super p-branes, 
i.e. K\ = 1. 

1) Single (d — l)-brane (N — l,n = 1) 

In the notations of sections 1 and 2 we have the solutions 



h{j 



Vrj 



,2C 



H 13 



A 



(0 



(3.16) 



where (3 = d/D. From these we find, G = 0, and for (dropping the r and A indices) 
becomes 

_ _ e (i) fj2a 2 -2+dd/2D 



So equation ( |3.9|) is satisfied, provided we have 



a 



1 - dd/2D 

which proves (|1.23|) again. Under this condition we have 



c {i) 



— e 



(0 



-AA' 







(3.17) 
(3.18) 

(3.19) 



2) Two orthogonal branes: di n d 2 = 8 
Using the solutions ( |1.29| ) we have 





(N = 2,n = 2) 








H 









I 
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^U(m) = Hi" e (At) e( m ) , A( M)(m ,-) = H 2 ~ e^e^') (3.20) 

where f3\ = d\/D and e^, ei m ) and e (m') are the alternating symbols corresponding to 
the volume forms of the subspaces: x^ , y™ and y™' respectively. All other components 
of the form-potentials are vanishing. Again we find G = 0, and the non- vanishing i^^'s 
become, 

■r \ = — e ^1 -"2 

p2(n)(m') _ _^^(^)^(rn')jj2a2 2 -2+d2d2/Djj2a 1 a2+5-d 1 d 2 /D ^ 2^ 



The only way for these to satisfy ( |3.9|) , is that all the exponents in (|3.21| ) vanish. This 
gives 

2 t ^1^1 
«i = 1 - 

a 2 2 = 1 - 




«i« 2 = 2 " « J (3.22) 

which are the same as the consistency conditions ( |1.21|) (with k\ = ±1). Under these 
conditions we see, in agreement with ( |3.9| ), that: = = const. The non- vanishing 

r(i) , 

c A s are 

c i(m)M = _ e ( M ) £ (m) 

c l M (m') = _ £ (M) e (m') (3>23) 

3) Type II A (2, 2) — branes at angles (N = 2, n = 1) 
The solution to the system 3 H 3 = 1 at two SU(2) angles in Type IIA theory, and its 
generalization for n number of such 2-branes, was originally conjectured by Myers et.al. in 
29| and later derived in |3T| . This is specially an interesting example for comparison with 



the results of this paper. Both of the above references have made use of an asymptoticly 
orthogonal Cartesian coordinate system (t, x±, x 2 , x 3 , x 4 ) for the world-volume directions. 
In these coordinates one of the 2-branes lies on the (xi, x 2 ) plane while the other is located 
on the plane (x'^a^), obtained from the first one by two SU{2) rotations with the same 
angle 9, 

I i ■ / id I i ■ \ 

x l + ix 3 = e (xi + txs) 

x' 2 + ix' 4 = e~ i9 (x 2 + ix 4 ) (3.24) 
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In these coordinates their solutions (in the Einstein frame) can be put into the forms 



ds 2 = E' 5/s [-dt 2 + (1 + ht + cos 2 9h 2 )(dx\ + dx 2 2 ) + (1 + sin 2 9h 2 ){dx 2 3 + dx\) 
+2cos9sin9h 2 (-dxidx 3 + dx 2 dxi)] + E 3/8 dz a dz a 
A = E~ 1 dtA[(cos 2 9hi + h 2 + sin 2 9hih 2 )dxiAdx2 — sin 2 9hi(l + h 2 )dx 3 Adx4 
+cos9sin9hi(dxiAdx 3 + dx 2 AdXi)] 

E l/A (3.25) 



where a = 1,...,5, hi and h 2 are independent harmonic functions vanishing at infinity, 
and E is defined by 

E = l + h + h 2 + sin 2 9h 1 h 2 (3.26) 

These solutions can be put into a symmetric form, relative to the exchange of the 
two branes, by transforming the coordinates (x\, x 2 , x$, x±) to the 'branes coordinates' 

. (^i , X 2 , X^ , X 

2) = (2/1,2/2,2/3,2/4), which even at z — > 00 do not constitute an orthogonal 
frame. This simple exercise gives 

d s 2 = sin- 2 9E- 5 / 8 [-sin 2 9dt 2 + (1 + h x + cos 2 9h 2 )[dy\ + dy\) + 

(1 + h 2 + cos 2 9h 1 )(dy 2 + dyj) - 2cos9(l +h ± + h 2 ){dy x dy z + dy 2 dy A )} + E 3/8 dz a dz a 
A = E^dtAlh^l + h 2 )dy\Ady 2 + [h 2 (l + hi)dy 3 Ady4 — cos9hih 2 {dy\Adyi + dy 3 Ady 2 )} 

e 2a„ = E l/8 ( 32 ^ 

where we have used the value of a = 1/4 for this case. From the above solutions we find 
G = — In sinO. We obtain the non- vanishing components of as 

^012 _ \ ^012 _ cos 2 9 

1 sin 2 9 2 sin 2 9 



COs2 ° ^034 = 1 

sin 2 9 1 sin 2 9 



014 _ cos9 014 _ cos9 

1 ~~ ' 9Q ' 2 ~~ • on 

sm z 9 sm z 9 

C0S9 „ncm cos9 



™ F 2 032 = =- (3.28) 

1 s%n 2 9 s%n 2 9 

which are seen to be constants, confirming our equation ( |3.9|) with c$ = 0. 

General formula for 

The main lesson we learn from the previous examples, is that the integration constants 
in the general ansatz ( 3.11]) are determined by relative orientations of the constituent 



branes. For this reason, we call them as the 'structure constants' of the brane system. 
A general formula for these constants, can be obtained using the boundary conditions. 
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For localized distributions (with d > 0), the metric and dilaton boundary conditions are 
written as 

/ir,-(z)| 2 ^oc = lij , C(z)\ z _ 00 = , ^(2)1^00 = ip = (3.29) 

where 7^ is a flat Minkowski metric with the signature (—,+,...,+) expressed in the 
(generally) non-orthogonal Cartesian coordinate system (x l ), and for convenience we have 
chosen (p Q = (c.f. section 1). The other set of boundary conditions are related to the 
form-potentials. These boundary conditions arise since as z — > 00, the gravity and dilaton 
fields become so weak as to decouple (to first order) from all the the form-fields. So at 
infinity, the set of form-fields propagate like an assembly of uncoupled fields (of gener- 
ally different degrees) within a flat background. These are described by the 'Maxwell's 
equations' : 

d*dA r = *J r (3.30) 
in the flat spacetime, with p r -form sources as 

J r = Y, /°a(«)ca (3.31) 

Here e\ and p\{z) denote respectively the asymptotic volume- form and the density func- 
tion associated with the X th distribution^ Solution to (|3.30|) is then: 

AU Z )= E txdMz) (3.32) 

d\ =Pr + l 

where h\(z) 7 s are harmonic functions defined by 

h x {z) = J d^ 2 z'p x (z')G d+2 (z,z') (3.33) 

Note that h\(z) — > as z —>■ 00, provided that p\(z) is a localized distribution. 
Viewing equation (|3.32j ) as the boundary conditions for the original problem, we find that 

^A^4(j)| 2 — 00 = &d x ,p r +i£\{i) (3.34) 
Then using the boundary values ( |3.29|) and ( |3.34|) into the equation (|3.12|) , we obtain the 



desired formula for c A 



r(i). 



$d x , Pr +iVi(i l3 ) e m (3-35) 

where 7 = |det(7y)|. An important caution for later applications of this formula is that 
it can be used only in a 'basis' of harmonic functions as is defined by (I 



11 This form of currents is read from the p r -branes er-model sources, coupled to the supergravity action. 
In the x l, s basis: e\ = e\^(dx l ) where the constants ev,-) for non-orthogonal x l, s may be other than 
0,±1. 
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Change of the harmonic functions basis 

For the equations describing A (//)'s (called as the 'embedding equations' ), we have an 
obvious symmetry under the permutations of H\s, as these equations do not involve any 
explicit dependence on //a's. This symmetry is a special case of a largest symmetry, i.e. 

H\ = a XK H K + b x , det(a XK ) ^ (3.36) 

which transforms the embedding equations covariantly (see (|3.56j )). This originates from 
the fact that any nonsingular linear transformation on a set of independent harmonic 
functions, transforms them to another such set. It turns out that //a's are similar to 
a set of Cartesian coordinates (called as the '//-basis' ) on a kind of 'geodesic' surface 
(called as the '//-surface' ) in the space of field variables |34j| . Symmetry of the embedding 
equations under fl3.36j ) is a reflection of the fact that, we have not any preferred origin 
and direction for choosing the //-basis on the //-surface. The transformation property of 
c x s under ( |3.36j ), can be read from their definition by fl3~P| ) to be 

■ r{l) - (a- 1 )^ (3.37) 



In section 4 we will find that the solution for a system of two p-branes at angles more 



easily described, in a certain basis which is related to the 'natural' one ( |3.33D by a suitable 
transformation (see ( f4.49| )). 

The reduced Lagrangian 

The gravitational and form-field parts of the total RL for (hij,C) and A r ^, are given by 
the equations ( |3.3|) and (|A.12|) respectively. Adding the dilaton term to them, the total 
RL will read 

£ = -e 2G [-Tr 2 (h- 1 9h) - ^Trih^dhf + 2(d + l)dC.Tr(hr l dh) 

n 

+2d(d+ l)(dC) 2 - (dtp) 2 - J2 e 2arlf (h ij )dA r {l) .dA r {j) } (3.38) 

r=l 

where d = d a , and G is defined as in Q3.4Q . Before discussing the equations of motion, 
some simplifications is in order: 

a) As in sections 1 and 2 we take the benefit of the constraint: G = Gq, which means 
that C must be replaced by 

C=-i(ln|h|-ln| 7 |) (3.39) 
where we have set Gq = l/21nypy( using the boundary conditions ( |3.29|) . Then by elim- 



inating C between ( |3.38j ) and ( |3.39| ), the simplified RL describing (hij,(p,A r ^) becomes 

n 

C = -l/4[J- 1 Tr 2 (h- 1 9h)+Tr(h- 1 9h)VV2(^) 2 -l/2^e 2a ^(^)M[ i) .M[ i) (3.40) 
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Demanding solutions to be on the 'surface' : G = Gq requires a consistency condition. 
This is obtained by a simple elimination of C from its equation of motion using ( 3.38|) 
and ( PE3gp . For d ^ this yields 

C = (3.41) 
which will be interpreted as the 'marginality' condition. 

b) As we have seen in sections 1 and 2, the number of actual metric variables for a given 
system, is often less than that of the hi/s. Denoting these variables together with the 
dilaton collectively by <p a (a — 1, M+l where <p M+1 = (p), the gravitational (including 
the dilaton) part of ( |3.40| ) will take a simple form as 

C G = 1/2^(0)^.9/ (3.42) 

Here £l a p's are known functions of a 's. This Lagrangian looks like a kind of (pseudo-) 



Riemannian metric on the space of a, s, with Vt a n as its metric tensor 34 



c) Let us consider, for a moment, H x s as the set of independent variables required to 
describe all the non-vanishing form-fields for the configuration of interest, i.e. Afy = 
A r ^(H x ). So the set of all independent variables is here ((f> a ,H x ), in terms of which the 
form-field part of (|3.40p is written as 

C F = ~Uxx'((f>,H)dH x .dH xl (3.43) 

with 

n 

U XX ,(<P,H) = ^^(h^dxAfodxAfa (3.44) 

r=l 

In the last expression the (known) (0 Q )-dependence of U X y comes through (</?, hij), while 

its (unknown) (i/ A )-dependence comes through ALy 

Using these points, the complete Lagrangian describing (0 a , H x ) becomes 

C = l/2n a p((f))d(f) a .d(j) /3 - l/2U xx ,{(j), H)dH x .dH x , (3.45) 

The equation of motion for H x from this RL is written as 

U xx ,d 2 H x , + dU \ X i . dH X i - l/2U m dH x ,.dH x „ = (3.46) 

where the symbol £ ,A' in the subscript represents the partial derivatives with respect to 
the explicit if^-dependences, wheres d x will represent the total i^A-derivatives. Now we 
enter the main three assumptions characterizing the set of marginal solutions as: 

1) d 2 H x = outside the region of branes distributions. 

2) Q 's can also be expressed as functions of H x s. 
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3) H x s constitute a set of independent harmonic functions. 
Using the above ingredients in the last equation, we obtain 



d x »U xx , + d x ,Uxx» - C/a'A'',a = (3.47) 

It is clear from ( |3.44[ ) that, the functional dependences of U xx >'s on (H x ) will not be 
known, unless we specify those of A^'s. We do not know these functions presently, but 
rather we have (from (|3.12|) ) the expressions for their derivatives d\ A r ^ as functions of 
(d a ). which using G — G^o, are written as, 

d x A\ i) =e- 2 ^c^\h t3 ) (3.48) 

Using this result in Q3.43j ), an explicit solution for Uw on the if -surf ace is find to be 

U X y(<P(H),H) = u X y(<f ) (H)) (3.49) 

where u xx >'s are given by 

n 

u xx ,((f>) = £e" 2 ^cf c^(M (3.50) 

r=l 

These functions are written as the sum of n homogeneous polynomials in the variables 
(hy), of degrees (p± + 1), (p n + 1), with the coefficients which are combinations of the 
structure constants. We will use, in section 4, these polynomial functions, to analyze the 
field equations of two p-branes at angles. 

In Appendix D it is shown that the above expression is in fact the solution of the equation 
( |3.47| ) for U XX ' on the //-surface. An important result of this analysis is that all the 
functions U X y (and so u xx >) on the //-surface can be derived from a kind of 'potential' 
function U(H), 

U xx >\h = u\\'\h = d x d x <U (3.51) 

which is unique up to any linear function of (H x ). (By this result, for solutions on the 
//-surface, we obtain: C F = — 1<9 2 [7. Then the £ = condition gives: Cq — \d 2 U). 

This property has an important role in the proof of integrability for a system of mul- 
tiply intersecting branes [pi"|. 



Embedding equations of the //-surface 

Having expressed every thing in terms of (</>"), it now suffices to set the equations deter- 
mining the 'embedding functions' <ft a = <f) a (H), to complete the solutions to the equations 
of motion . These may be found simply by starting from the Q 's equation of motion 
derived from ( |3.45| ): 



n aP (<P)d 2 <P + Y a p^(j>)d^.dcjP + 1/2C/aa' iQ (0, H)dH x .dH xl = (3.52) 
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where r a/ g 7 is defined in terms of tt a p as 

= 2(^/8,7 — ^/?7,a + ^a,p) (3.53) 

Of course the set of equations ( |3.52| ) can not be solved unless we have specified the explicit 
forms of Uw's as functions of (4> a ,H x ), which is not the case at this point. But we can 
use a surprizing property of these functions, 

[/ AA> (0,iJ) = -« AA ,, a (0) (3.54) 

(which gives the derivatives U XX i )Ce as given functions of (<f) a ) on the if -surface, see Ap- 
pendix D), to rewrite ( |3.52j ) as 

tt a p(<P)d 2 ^ + r afh (<j>)d<fJ , .d<F> - l/2u X y ta (<p)dH x .dH y = (3.55) 

which is now a well-defined equation. Looking for solutions on the if-surface, is equivalent 
to using the above three assumptions in this equation. This leads to 

n a p(<j>)d x d x ,<j>P + V aPl {(j>)d x ^d xl( jP - l/2u xx>> M = (3.56) 

This equation for the 'embedding' of the ii-surface, is similar to a a forced geodesic 
equation corresponding to the metric defined by 0^(0). As stated earlier, this equation 
does not involve any explicit iiA-dependences, and is symmetric under the transformation 
Q3.36j ) provided u xx ^s transform as 

uw = (a _1 ) K A(a _1 ) K 'A'M KK ' (3.57) 

One can check, using the definition of maa', that this is in fact the case. 
It is instructive to note that equations ( |3.56| ) can be put into the form of a true geodesic 
equation (i.e. without force), if one extends the configuration space to <p A = {(p a ,H x ) and 
defines its metric by 

dS 2 = n AB d(f) A d(f) B = Vt^d^dcj) 13 - U xx ,dH x dH x , (3.58) 

Then, of course, the if -surface will be a null-geodesic surface in the space of A 's. 



Constraints 

The formulation described above provides a general framework for studying the structure 
and properties of the marginal solutions. In particular it has the potential to give the 
constraints of sections 1 & 2, and their generalizations to the systems of branes at angles, 
in a systematic and exclusive way. We postpone such issues to a forthcoming article |Hj . 
Here we give heuristic derivations for the generalized versions of the constraints ( |1.10| ), 
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for arbitrary systems of (multiply-intersecting) branes at angles. The justification for 
use of these constraints is that they have simultaneous solutions with the equations of 
motion . This can be verified partially, using the solutions in sections 1, 2 and 4 of this 
paper. Generally the constraints for a marginal configuration, consist of the two following 
categories: 

a) Extremality condition: 

This constraint is somewhat related to the 'marginality' condition, which states that 
the binding potential energy of such a system, for arbitrary positions of its constituents, 
vanishes. In the language of fields, this means that the Hamiltonian density of interactions 
7i must vanishes everywhere. So on a time- independent marginal solution with 7i = —C, 
we must have C = 0. Noting the C's equation of motion as 

2(d + l)d 2 e 2G - l/2dC = (3.59) 

we see that a marginal solution satisfies 

d 2 e 2G = (3.60) 

i.e. the combination e 2G of the metric components must be a harmonic function. The 
most natural choice is to take G = Go, or equivalently 

In J\h\ + dC = In (3.61) 

It must be remembered that the C = is not a single equation in this formulation; it 
contains a set of component equations. This is because, by evaluating the RL on the 
iJ-surface we obtain 

C = l/2[n aP {4>)d x 4> a dw<t> P - u xx ,(<P)]dH x .dHy (3.62) 

By independence of H x s, the C = equation implies that 

£xx> = n a Mdxrdy^ - «aa'(0) = (3.63) 

which introduces a set of ^N(N + 1) equations (of I s * order), satisfied by the (M + 1) 
functions <f) a (H). If we have enough number of constraints, we can combine them with 
these equations to solve for a (if)'s; a method which will be used in the next section. 

b) No — force conditions: 

There are another set of constraints, though not seemingly expressed as the no-force 
conditions, but are equivalent to the latter. The no-force conditions corresponding to a 
configuration of iV distributions, are written as a set of N constraints relating the dilaton, 
metric and form- fields components together (see section 5). By expressing (formally) the 
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form-fields in terms of the harmonic functions, we can re-express these conditions as 
relations between (0 a ,if^)'s. We obtain such relations in the following using a heuristic 
way of derivation: We consider a special case where all the charge distributions except 
the X th one are vanishingly small, i.e. p K (z) — > for every k ^ A. The A r, s solution for 
such a single distribution is known to be 

A r = -d^^Hx-hx (3.64) 

where Hx = 1 + h\ with h\ as in Q3.33j ). The corresponding part of the RL is calculated 
to be 

jC fx = l/2e- 2 ™\dH x - 1 ) 2 (3.65) 
where the function F x (<f)) is implicitly defined by 

It is easy to check that the Hx's equation of motion from ( 3.65Q becomes a Laplace 
equation, only if we have 

F A (0) = -\nH x + const. (3.67) 

This gives the desired constraints, which are claimed to be valid for a general situa- 
tion. Note that in the special cases where af's are chosen to be the branes' world- 
volume coordinates, exu) becomes the alternating symbol corresponding to the coordi- 
nates on (d\ — l)-branes' world-volume, and the definition (|3.66| ) simplifies to e 2Fx< ^> = 
e W+°w)\det(h*1)\ i)jedx . 



4 Application to a system of branes at angles 

An illustrative example of the formulation of section 3, is provided by considering the 
distributed system: d\ D c?2 = 5 at non-trivial angles. In the branes' world-volumes 
coordinates (as defined in section 1) the ansatz for the metric is written as 



ds 2 



e 2B (z) dx n dx ^ + e 2Bi(z) dy m dy m + ^{z) dy m> dy m> 

+2e Bl{z)+B2{z) q{z) lmm ,dy™dy™' + e 2C{z) dz a dz a 



(4.1) 



This differs from its orthogonal analogue, equation (|1.2j ), by an off-diagonal term, which 
appears due to the non-trivial angles between the branes. In this expression q is some 
new degree of freedom and 7 mm ''s are S1S2 parameters defining the angles. The internal 
(world-volume) part of the metric tensor here is written as 

fe 2B ^ u 









V 



qe Bi+B * lnm , 



qc Bl+B2 -lmn> 



(4.2) 
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Rotating the internal coordinates on the individual branes changes the matrix 7 = (7 mTO /) 
by a transformation of the form: 7 — > A7_B T , where A G and -B G SO^)- 

However, such different 7's describe geometrically equivalent angular situations. The 
SO(5i) x SO(52)-rotationally invariant parameters are the angles: {9i, 9 m }, defined 
through the diagonalization of the matrix 77 71 (see Appendix C), whose number is 

m = Min(6i, 5 2 ) = Min{d x , d 2 ) - 5 (4.3) 

This is smaller than the total number of r y mm ' , s. 

The expression for £ G 
With the aid of the formulas A. 12 to A. 20 of the Appendix, one obtains the RL, 



£g = f(q) 1/2 e 2G ^{l/2n a/3 (q)dr.d^ + a{q)a a d<j> a .dq + l/2b{q){dq) 2 } (4.4) 

where (fi a = (B ,Bi,B 2 ,C,(p) and G(<fi) is defined by Q1.6Q and (|1.7|) . Further we have 
defined the functions of q as 

f(q) = det(l - g 2 77 T ) 
a(q) = f(q)/2f(q) 
6(g) ee f"(q)/2f(q) 

& a p{q) = n aj3 + u a upqa{q) (4.5) 

with a prime denoting d/dq. The constants (u a ,a a ) in these formulas are defined 

u a = (0 1 -1 0) 

a Q = (25 2<Ji-l 2<J a — 1 2(d+l) 0) (4.6) 



The matrix Q a/ g in [4.5| is defined by ( |1 . 7| ) . According to the above expression, the explicit 



angular dependence of Lq comes totally through an even polynomial f(q) of degree 2m 
, whose roots {q r } characterize the angles between the two branes via qf, = l/cos 2 8 r . In 
the orthogonal case with •y mm ' = 0, we have /(g) = 1 and consequently the expression for 
Cq returns to that of eq. (|1.4|) 



Intersection rules revisited 

A natural question is that whether the intersection rules of the orthogonal branes remain 
to be valid even for branes at non-trivial angles ? We argue that the main equation de- 
termining these rules can not be angle-dependent and so the answer to this question is 
positive. First of all, since the mechanism by which we obtain a marginal solution in this 
paper, is essentially the same for the two cases, we expect that any non-orthogonal solu- 



tion of this type to be connected to the orthogonal solutions ( 1.29 ) continously through 
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the parameters (9i,...,9 m ). So in particular, we should expect to require three consis- 
tency conditions analogous to the equations (|1.21| ), relating the couplings, dimensions, 



and possibly angles. But the last possibility can not be true, as if the counterparts of 
equations (|1.21|) include angles, then either the expressions for (0:1,0:2) or the relation 
arising between dimensions (counterpart of ( |1 -28j )) must contain angles. Neither of these 
possibilities are physical. Firstly (oi, a 2 ) can not be angle- dependent, as they are param- 
eters of the model itself which can not depend on the geometry. On the other hand a 
relation between dimensions and angles implies a quantization of angles which has not 
any physical origin. Therefore in what follows we will assume that the equation ( |1.28|) for 
super-p-branes is always applicable. 

Derivation of the solutions for (p,p)-branes at angles 

To present an application of all the previous general formulas, we show here the stages 
of derivation for the solutions of two 'identical' branes at angles. So in all the previous 
formulas, we put di = d 2 = d and a\ = a 2 = Oo- 111 this case we have only a single 
form-field A of degree d Q . Applying ( |1.28|) to this system, we obtain 



5x = 5 2 = 2 
S = d - 2 

d = d Q + 2 (4.7) 

The relative orientation in this system, (generally) is described in terms of two indepen- 
dent angles, say (6,8'). A simplification occurring here is that, both the subspaces of: 
(2/™) and (y™') become 2-dimensional. So 7 mm ' becomes a 2 x 2 matrix, which after 
diagonalization takes the form 

(cos9 \ 
7mm' = _ „,} (4.8) 

\ cosO J 

Hence, putting: y™ = (2/1,2/2) and y™' = (2/3,2/4), the metric (fOJ) reduces to 



ds 2 = e 2Bo ^dx»dx, + e 2B ^)(dyl + dy 2 2 ) + e 2B ^ z \dy 2 3 + dyl) 

+2e Bl{z)+B2{z) q{z) (cos6d yi dy 3 + cos6'dy 2 dy 4 ) + e 2C{z) dz a dz a (4.9) 

We assume in the following the asymptotic flat boundary conditions (provided d < D—4), 

a Uco = , gUco=l (4.10) 

The i/-basis and the structure constants 

We begin our analysis by computing the structure constants in the basis h\, defined 
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by the equations ( |3.35| ) and ( |3.33| ) respectively. This requires the asymptotic form of hij, 
which by (|4.9| ) and ([4.10 ) is written as 



7ij 





















1 





cos9 











1 





cosd' 





cos9 





1 





V o 





cose' 





1 



(4.11) 



where the indices: r,s = 1, ...,4 refer to the coordinates (yi, 3/4). Then the non- 



vanishing structure constants cf^ rs according to ( |3.35| ) and ( |4.11|) are written as 



6 C A 



(4.12) 



where 7 = dct(^ rs ), and and eA( rs ) (with A = 1, 2 and r,s = 1, 4) are respectively 
the (p — 1)— and 2-dimensional alternating symbols such that: e '"^ -2 -' = —1 and £1(12) = 
e 2(34) = +1. A straightforward calculation using fl4.11| ) and Q4.12j ) then gives 



-12 



,3 1 



„32 



,11 



,-,31 



-.12 



r 14 

^32 



A 

\cos6cos6' 

\cos9 

\cos6' 



(4.13) 



where A 



l/sin8sin6' and all other c^ s 's are vanishing. Note that for 9 = 9' 
, as is so for the (2,2)-brane solution ( |3.27| ), the resultant values of c™ exactly reproduce 
the constant values of , which were given by (|3.28| ) . We now introduce a new if -basis 
by 

H x = l + a XK h K (4.14) 

with a transformation matrix (a\ K ) as is determined below. The constants c r x s in the new 
basis (H\) then is related to the old one (h\) according to ( |3.37| ) by 

c\ s = (a- 1 )^ (4.15) 

To determine the coefficients (a K x) we demand that: 

i) «kA = cl\k so that our solution to be symmetric in (Hi, H 2 ), as is so for (hi, hz). 

ii) Cg 2 = cf 4 = 0, which is necessary for the integrability of our equations. 

iii) For 9 = 9' = n/2 we must have H\ = 1 + h\. 
These together fix (a K \) as 

/ 1 cos9cos9' \ 
, cos9cos9' 1 / 



(4.16) 
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Using this matrix in ( f4.15| ), we obtain 

c{ 2 = c^ 4 = A 
cf = t\ 4 = \v 

c} 4 = cf = Xu' (4.17) 

where v = cos0sin 2 0' / (1 — cos 2 0cos 2 0'), v' = cos0' sin 2 0/ (1 — cos 2 6cos 2 9 ) and as in the 
past A = \/sin0sin0'. 

Explicit expressions for uw^s 
Knowing the values of the structure constants in the H\ basis, we are able to obtain the 
functions u\y (0, q) , also in this basis, using the formula ( |3.42| ) which reduces to 

uaa' = -e- 2id - SBo+a ^w X y (4.18) 

Here a hat is used to distinguish between the analogous quantities in the orthogonal and 
at angles cases (see (|4.26|) below) , and wxy's are homogeneous polynomials of degree 2 



in (h rs ) defined as 

w X \> = c r x s c r y s 'h rr ,h ss < (4.19) 
Using the ansatz ( |4.9| ) for the metric and the values of the structure constants as in ( |4.17|) 



we can expand the above formula to obtain 



wu = A 2 (/ii 1 + ahnh 3 3 + 2phuh Vi + rh 2 



13J 

w 2 2 = A 2 (/43 + <7h u h 33 + 2ph 33 h 13 + rh 2 13 ) 

w X2 = \ 2 {^h n h 33 + p'h u h 13 + 2p'h 33 h 13 + v{l + a)h 2 13 ] (4.20) 

The constants appearing in these formulas are defined in terms of the angles (8, 9') as 

cosO' 



v = 



cosO 

p = vv + v' , p = vv' + V 

a = v 2 + v' 2 , t = Ivvv' (4.21) 



where v and v' have been defined in terms of {0,0') below (|4.17|) . We now introduce the 
'homogeneous' variables (g, s), 



Vh n h 33 V h u 



which are invariant under the re-scalings of (h rs ). Then using ( [4.18|) and ( |4.20| ) we obtain 
an expression for u\y as 



u\\> 



_ A 2 e2 (-W +BA+5v -M /AA , (giS) = -X 2 e P ^'fMg,s)/f(q) (4.23) 
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where f\y(q, s)'s are three polynomials in (g, s) defined as 



f u (q,s) = (1 + as 2 ) + 2psq + rs 2 g 2 

f 22 (q,s) = (l + ^) + 2^ + V 
s s s 

f l2 (q,s) = - + p'( s + -)q + v(l + ( j)q 2 (4.24) 

V s 

In (|4.23|) , G and F\ are generalized versions of G and for the case at angles, i.e. 

G(<j>,q) = G(0) + l/21n V / 7(g) 

F x (<f>,q) = F A (0)+l/2m//(g] (4.25) 

The explicit expressions of the functions appearing here are 

G{4>) = l/2(5B + 2B 1 + 2B 2 + dC) 
Fi(0) = 1/2(5 B + 2B 1 - 2B 2 — dC — 2a p) 
F 2 ((j)) = l/2(5B -2B 1 + 2B 2 -dC -2a <p) 

f(q) = (l-g 2 )(l-t;V) (4.26) 
The orthogonal situation is recovered by putting: /(g) = 1 , f\\>(q, s) = 5\y, and A = 1. 



Constraints 

The two classes of the constraints ( |3.61| ) and ( ^.67[ ) for this case take the forms 



G(0,g) = l/21nV7 

F x (<p, q) = - In H x + 1/2 In (4.27) 

where the constants on the right of these equations are obtained from the boundary 
conditions ( f4.10| ) and that H\—>-l as z^oo. These constraints relate the six variables 



i a ,g) and the two H\S. Equations ( |4.27|) , compared to their orthogonal analogues 



in (II.IOQ , have the additional term 1/2 In y /(g) /*y. This suggests that one may find a 
consistent simultaneous solution of these and the equations of motion , by taking 



<j ) a = -C x X x + V a \n^f(q)/ 1 (4.28) 

where X\ = — lnH\ and (£,\,r) a ) , s are a set of undetermined parameters (compare to 
( |1.17p ). Replacing this in ( |4.27| ), and taking (Xi, X 2 , In a/7) as linearly independent func- 



tions, one obtains nine algebraic equations for the parameters. This consists of the equa- 
tions ( 1.18|) and the following equations 



g.y/ = f x . V = -1/2 (4.29) 
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where g and fA are defined as in ( |1.7| ). To completely determine these parameters, the 
consistency with the equations of motion must also be required. However, by considering 
the 'single distribution' limits one can verify that the values of £°'s must be independent 
of the angles and in fact they do depend only on the dimensions of the respective p-branes 
. This means that £"'s are already determined as in (|1.20|) (with k\ — k| = 1). The values 



of r] a, s are then determined using the £ = condition (see below). 

The £\y = equations 

As was stated earlier, all the field variables: (q, (fi a , A^ rs ) must be expressible as functions 
of H\S. These are not independent functions; since by ( |3.12j ) and ( [4 . 2 8| ) , a knowledge of 
q{H) determines all other functions of H x s. To determine the former we use the £ = 
condition, which by (|3.63|) gives 



e 2G [£l a pd x <t> a d y( j>P + a a a(d x <p a dx>q + dy^dxq) + bd x qd x ,q] + e F ^'f X y/f = (4.30) 
where dx = d/dHx and 

s = e ~(B,-B 2 ) = e -u, a r ( 431 ) 

Using (|4.27| ) and replacing H x with X x , eq. ( f4.30|) becomes, 



[aHrii'hl + 2a.Ti)+b\d x q/^ 

(4.32) 

where dx = d/dXx- This can be put into a more convenient form, 

(d x q - a x (q))(d x ,q - a x >(q)) = b xx ,(q, s) (4.33) 

where we have defined; 

C 

bxx> = --(t>fex> + fxx>/f) 
c 

c = b + (ririflri + 2<j.t/) (4.34) 

Note that in ( [4.34|) /, a, b, c, ax, ^2 are only rational functions of q, while f xx > and bxy are 
rational functions of q and polynomial functions (at most of degree 2) of s and 1 fs. By 
subtracting the A = 1, 2 equations in |4.27| , we obtain 



Hl /H 2 = e - 1 /2(x 1 -x 2 ) (4>35) 
The integrability conditions 

The three partial differential equations ( f4.33| ) (for A, A' = 1, 2 ) have a solution if, 
i) a consistency condition and zz)an integrability condition are satisfied. 
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To see this, consider first the two equations with A = A' . These determine d\q's as 
functions of (q,X x ), 

d\q = a x + ex\fb\\ (4.36) 
where E\ = ±1, upon which eq. ( |4.33| ) requires the consistency condition: 

bnb 22 = (b 12 f (4-37) 

The integrability condition follows from: d\(d 2 q) = d 2 (d\q), which using ( |4.36| ) gives 



exJbn (4a'b 22 - 2ab 22 + sb 22 ) - e 2 \^b 22 {\ab xx - 2ab' n + sb n ) = 2(b n b' 22 - b 22 b' u ) (4.38) 



where prime means d/dq and dot means d/ds. By symmetry we have set ai = a 2 
Comparing (|4.38|) to an equation in the domain of Z of the form0 



a. 



ay/x + b^fy 



(4.39) 



one concludes that b\\(q, s) must be written as the square of a polynomial p\(q, s) in s or 

1/8, 



bxx(q,s) = [px(q,s)] 2 (4.40) 

This statement, as will be seen, leads to important implications about the angles and other 
parameters of the solution. Although logically fl4.40|) is only a necessary integrability 



condition for ( |4.36| ), in the case at hand it provides the sufficient conditions as well. 
Explicit expressions of &aa''s follow from their definitions by (|4.34 ) restricted to the case: 
d PI d = do — 2 via the equations (|4.24| ) and ( [4.26|) . In this calculation we encounter 
various combinations of the parameters {£,\,n a ), all of which except one are found to be 
calculable from the equations ( |1.18|) , ( |1.19|) and ( f4.29| ). For later reference we list them 
below 



-S 



XX' 



V^x = "I 

= 1/2 



e A «e A , = 1/4(2^ - 1) 

rfu£ x = rjurj = 
r/ fir/ + 2a. Tj = k 



(4.41) 



Only the last combination (k) remains unknown which must be determined using the 
condition (|4.40|) . 

The resultant expressions for bu and b 22 become 



hi(q, s) = b (q) + h(q)s + b 2 (q)s 2 
b 22 (q, s) = b (q) + 61(g)- + b 2 (q) — 



(4.42) 



12 Mathematically this comparison is possible, since the set of polynomials (in (s, 1/s)) , as that of the 
integers, with ordinary addition and multiplication, constitute a 'ring'. 
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with g-dependent coefficients as 



a 2 a 1,„ 1, 
60 = 4?-^ + c (1 "7^ 



2pq 



h = (4.43) 



fc 
-f 
/c 

Now for b\\(q, s) to be the square of a polynomial (in s or 1/s), we must have 

bl(q) = Ab (q)b 2 (q) (4.44) 
For this to be held identically the relations among coefficients must be 

v 2 = l 

a = r = p 2 /2 

k = — 2 (4.45) 

which leads tc0 

cose' = cos9 (4.46) 

■qn-q + 2a.r] = -2 (4.47) 

The intersection rule ( |4.7|) , combined with relation (|4.46|) , leads to the interesting conclu- 
sion: 

A pair of identical supersymmetric branes in any dimension make a a marginal configu- 
ration, if the two branes intersect at two equal or equivalently at SU(2) angles. 
The existence of a configuration of two intersecting branes at SU(2) angles is a result of its 



unbroken 1/4 supersymmetry first established in [28]. A similar conclusion was reached 
in the context of string theory calculations in ||33|| , where it was shown that the inter- 
action potential between two D-branes at SU(2) angles identically vanishes. Our result 
establishes their conclusions in a purely field theoretic context. The advantage of such a 
derivation is its independence of the high-energy model of the theory and in particular of 
the spacetime dimension. 

We will return to the other result, equation ( |4.47|) shortly. It is suitable here to emphasize 
on another feature of the integrability conditions which is only implicit in our calculations, 
i.e. their role in fixing the suitable if-basis for using with the constraints ( [4.27[ ). This 
relation becomes clear, recalling the role of u>aa''s (see ( [4.18 )) in our equations. These are 



13 We must discard the possibility of cosO 1 — —cosQ as will be seen below. 
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2 nd degree homogeneous polynomials in h rs variables, and have six terms generally. How- 
ever, by choosing a suitable if-basis, we can remove two of these terms simultaneously. 
For example the expression for wu has the form 

w n = (cf) 2 (h 33 ) 2 + 2c\\c\ 4 + cf )h 33 h 13 + ... (4.48) 

where the doted terms do not contain any powers of cf 4 . We can suppress the first two 
terms by choosing an if-basis for which cf 4 = 0. As a result the terms of higher degree than 
s 2 (i.e. qs 3 and s 4 ) will be removed from fn(q, s), so that finally a quadratic expression 
(in s) for &n(g, s) is obtained. If we choose a basis for which cj, 2 = 0, a similar result 
for 622(9? s ) emerges. These two conditions (among the others stated below ( |4. 15| ) ) then 
determine, as we have seen, the matrix (a K \) defining the suitable basis. The important 
point here is that without taking: cf 4 = c}, 2 = 0, our equations Qi.36f) for q(X) would not 
be integrable. The suitable basis, using ( (4.16| ) and 9' = 9, is therefore 



Hi = 1 + hi + cos 2 9h 2 

H 2 = l + cos 2 9h 1 + h 2 (4.49) 

The solutions for q(X) and 4> a (X) 
Using 9 = 9' in our equations causes considerable simplifications. A straightforward 
calculation in this case, using fl4.24p and ( [4.34| ), shows that 

/(?) = (l-<? 2 ) 2 
2q 



1 — q z 

ai(q) = a 2(q) = - 



g(l + 9 2 ) 
2(1 ~q 2 ) 



&22(g,s) 



(as + q + asq 2 ) 2 

(1 -q 2 ) 2 
(f + q + - s q 2 f 



612(9,5) = r, 2^ — ( 4 - 50 ) 



(1 -q 2 ) 2 
(as + q + asg 2 )(^ + q+ f q 2 ' 

iw 

where a = cos9 / (l+cos 2 9), and we have used the numerical relations (|4.41|) with k = —2. 
It is clear from these expressions that the consistency condition (|4.37|) is automatically 
satisfied. Note, however, that if we chose the other possibility from (|4.45|) : cos9' = —cos9, 
then the consistency condition would not be satisfied. Using these in the equation(s) ( |4.36| ) 
leads to four possibilities for solutions of (diq,d 2 q), from which only the E\ = e 2 = — 1 
case is integrable. This corresponds to 

A 1 
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d 2 q 



i a 



(4.51) 



where s = e 1 / 2 ( Xl x ^ as in ( |4.35|) . This system with the boundary condition: q(0, 0) 
cosd has a unique solution, 



q = a[2cosh(X 1 ~X 2 )/2-Pe {Xl+X2)/2 ] 



a(H x + H 2 -P)/JH x H 2 



(4.52) 



where (3 = sin 2 8. 

The solutions for Q (X)'s, using this solution for q(X), are already expressed as in (|4.28|) 
up to the unknown parameters n a . To determine n a, s (as well as £"'s in principle), we 
have to set the relations between parameters, required by the consistency of the equations 
governing ((p a (X) , q(X)) and the ansatz ( 4.28| ) for (f) a (X, q). Four of these relations have 
been given previously by ft4.29|) and (|4.47| ). For obtaining a fifth one, we use the embedding 
equation of the dilaton, 



(d x dx> + 5 X \>d x )p + aofxx'/f = 



(4.53) 



where have used the constraints ( |4.27| ) to eliminate the exponentials. On the other hand, 
the dilaton (i.e. a = 5) component of ( |4.28 ) is written as 



-a (X 1 + X 2 )+r ] 5 \nJf(q)/ 1 



(4.54) 



where we have put £f = ^| = a using ( |1.20| ). Inserting this ansatz in ( [4.53| ), a new set 
of equations involving only q(X) results. For these to become identities, equations ( |4.51| ) 
require that: rf = a . Putting this value in the equations (|4.29|) and (|4.47| ) and solving 
for n a, s, then gives 



do _ do _ do do 



(4.55) 



This completes our solution (|4.28| ) for (p a (X). 



The solutions for A{X) 
For purely electric-type branes only the 'world- volume' components of the form-potentials 
survive (see (|3.1|)), 

A = A^ysidx^Ady'Ady 3 (4.56) 
Applying the formula ( 3.12| ) to this case, we obtain 



d x A 



'A /1 (/^)rs — e (/i) c A ll>rr' n s 



T ' s 'h ,h rp 2(SB -G-a <p) 



(4.57) 
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With the help of ( |4.27j ), and in the form- notation, eq. Q4.57 ) is written as 



dA {fl ) rs = e^)U rs (4.58) 
where cj's are a set of 1-forms (on the if -surface) defined by 

uirs = ^H x - 2 dH x (4.59) 

where /x( rs )'s are functions of (q, s) with the non- vanishing members: 

2 

q q 

/i(i2) = 1 + 2asq / 2 (i 2 ) = 2a— + — 

,3„ , „2 2 J- 1 , r, _ 5 



/i(34) = 2as g + sg / 2 ( 3 4) = l + 2a 



/i(32) = /i(i4)=as 2 + sq + as 2 q 2 / 2 (32)=/2(i4) = ^ + ^ + ( 4 -60) 

For the equations (|4.58|) (for A^w's) t° be integrable, the forms kv s must be exact, 
which can be checked by expressing them totally in (q, s) variables. Therefore u's can be 
integrated to yield, 

e rs = - / uj rs (4.61) 



The non- vanishing e rs 's are; 



ei2 



e 3 4 



1 + 4 + ^ 

s z as 

(3(1 - g 2 ) 
1 + s 2 + f 



e 32 = eu,^ (4-62) 
Finally the solutions for the components of ^4 are written as 

A^rs) = -£(^e rs (4.63) 

Summary of the solutions 

In the following we present the final forms of the solutions in terms of the variables 
(H\, q), without inserting the if-dependences of q in them. This choice makes them more 
transparent. First, the expressions for ds 2 and are found from ( |4.28| ), using ( |1.20| ) and 
(PI), as 

/ / \ -do/B 

ds 2 = (y/(<z)/7tfi#2j [dx^dx. + H^dyl + dyD+H^dyj + dyl) 

+2q y /H 1 H 2 (dy 1 dy 3 + dy 2 dy 4 )} + [y] f{q)hH x H^ ,D dz a dz a (4.64) 
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Ola 



lf{q)hH x H 2 

Second, the expression for A is found, from (|4.62| ) and ( 4.63]) , as 
I 1 , i , o 1 . s , a q 



(4.65) 



A 



f(q)h 



{dx»)A[^-dy l Ady 2 + ^-dy*Ady i + 



(dy 1 Ady* + dy 6 Ady 2 )} (4.66) 



In these formulas, if -dependence of q is given by ( |1.52j ) and f(q) is defined in ( |4.50|) . The 
harmonic functions (Hi(z), H 2 (z)) are related to the charge densities (pi(z) , p 2 (z)) via 
( [3.331 ) an d ([4.49|) which altogether may be written as 



H x {z) 
H 2 (z) 



1 + J d d+2 z'G i+2 {z, z')[ Pl {z') + cos 2 6p 2 (z')] 
1 + [ d* +2 z'G d+2 (z, z')[cos 2 e Pl (z') + p 2 (z')\ 



(4.67) 



By expressing (H^H^ in terms of (hi, h 2 ), it is easy to check that for Type IIA (2,2)- 
branes with D — 10 and do = 3, this solution exactly matches that of p9f expressed 
as ( p. 27] ) previously. In particular the function E in that solution is nothing but: E = 
J f(q)/"fHiH 2 . The general case with arbitrary D and p had not been reported earlier. 



5 No-force conditions 

No-force conditions are a set of constraints arising naturally when we consider the marginally 
stable (static or stationary) configurations of p-branes . We have seen how a class of field 
theory solutions for these systems are constructed based on a very special set of con- 
straints corresponding to the extremality and no-force conditions. The interpretation of 
the former was given in section 3. In this section we give a precise interpretation of the 
latter, which justifies its name. We start from a formulation of the no-force conditions 
similar to that of f25fl . In principle to find such conditions, one has to single out every 



constituent (d\ — l)-brane from rest of the branes, and demand its 'equilibrium' condi- 
tions under interactions with the others. Equivalently we may look for the equilibrium 
conditions of a (d\ — l)-brane 'probe' , with nearly zero mass and charge, situated parallel 
to the similar distribution within the system. Since the spacetime geometry along the 
world-volume of the brane in equilibrium is homogeneous, we can decompose the metric 
as 

ds 2 = h iiV (y)dx il dx v + k mn (y)dy m dy n (5.1) 

where (x M ) and (y m ) stand for the coordinates parallel and transverse to the world-volume 
of the (d\ — l)-brane respectively. The action governing the dynamics of the brane is the 
DBI action. Since our branes have no boundaries, we can truncate DBI actions such that 
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all the internal gauge fields of the branes vanish |H2|. Therefore the static gauge action 
of a (d\ — l)-brane with X M = £ M , Y m = Y m (£) is written as 

S X [Y] = —T\ J d d ^{e- aMY) [-det (h^(Y) + k mn (Y) d fJ Y m d v Y n )] 1 ^ + k^s^A^Y)} 

(5.2) 

Here T\ is the brane tension, e x ^ the c^-dimensional alternating symbol, and ^4a(^) rep- 
resents the pull-back of the c^-form potential on the world- volume of the (d\ — l)-brane. 
In (|5.2|) we have introduced an unusual factor K\, which must be the same as the mass to 



charge ratio ( 1.24 ) of the brane (see below). For = ±1 we recover the usual action of 
super p-branes . Expanding this action in powers of the velocities: (d^Y m ) and keeping 
only the leading order terms, we obtain the static potential between a (d\ — l)-brane 
probe and the complete brane system as 

V X (Y) = T\ (e- a ^y/h x (Y) + ka-^AmOO) (5-3) 

where h\ = \det{h,j, v )\. In the equilibrium state, we must have V\{Y) = constant. 
Eliminating this constant by a suitable gauge transformation on A\, the no- force condition 
for a (d\ — l)-brane is written as V\ = , alternatively 

e -a xVh l/2 = _ Kx -l £x M Am (54) 

Note that the above equation gives explicitly the component of a c?A-form potential parallel 
to (d\ — l)-brane, as a function of the dilaton and the (determinant of the) corresponding 
metric components. This can be compared to an analogous but different result, which 
was given in the form of equation ( |3.12| ) previously. As the the number of the equations in 
( |5.4|) is equal to the number of constraints (p.67|) , N, we may tend to guess that these two 
sets of equations are equivalent. This is in fact true for all the configurations considered 
earlier, as will be seen below. 

The case of N orthogonal branes 

In this case by (|2.2|) , for all the branes within the X th distribution, we have 

A Mp) = e x(p) e x * 
e -w h y* = expC£hiBi-a x <p) (5.5) 

i 

So noting e Xx = KxHx -1 , equation ( |5.4| ) gives 

^6xiBi — a\<p=-1jiHx (5.6) 

i 

Note that cancellation of the factor n\ in this equation, occurs only if one introduces the 
factor ka _1 in the action ( |5.2| ). This confirms the result of section 1 that k\ is the mass 
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to charge ratio of a (d\ — l)-brane. We can put (|5.6| ) into a more familiar form, using the 



extremality constraint (equation (]3.61|) ) 



(5.7) 



Combining the last two equations gives 



Fx{<l>) 



l/2dC-a X ip= -In if* 



(5i 



This is exactly the second constraint in ( 1.10|) or its generalization to (|3.67j ). Clearly this 
derivation does not rely on supersymmetry properties of the theory, as it does not rely on 
choosing k x = ±1. In this sense, the marginality of a solution is not necessarily a result 
of its supersymmetries. 



The case of (p, p)-branes at angles 

For comparing with section 4, we consider only the case of super p-branes. In this case, 



using the ansatz (f4.2p for hij and the solution 



for A, we obtain 



A x = H x - 1 [f(q)/ 1 ]- 1/2 



where we have introduced A x = e\^> Auy So 



1/2 _ 2B X +SB - 



a tp 



takes the form 



2B X + 6B - a Q <p = - \nH x - In Jf(q)/~f 



(5.9) 



(5.10) 



Again the extremality constraint (first equation in Q4.27 )), combined with this equation 



gives 



F A (0) + 1/2 In yjf(q) = —h\H\+ 1/2 In 7 
This represents the second constraint in (|4.27|) , as was expected. 



(5.11) 



Long range potentials 

The (short and long range) potentials between D-branes in the spacetime of dimension- 
ality D = 10, due to exchange of (massive and massless) closed superstring states, have 

|, B7|. However there has not been a general pre- 



been calculated in many places |35 



scription for such calculations in spacetimes of arbitrary dimension, as a quantum theory 
of strings in other dimensions does not make sense. Nevertheless by studying the long 
distance behaviour of the 'effective' potential ( |5.3|) , we can specify the contributions from 
the various (massless) bosonic fields in our (low energy) model to the total long-range po- 
tential between a brane probe and a system of distributed branes. This will give physical 
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interpretations to the various consistency conditions, arising in the discussion of marginal 
solutions (e.g. eqns. flOT )). To this end, we first note that the various dynamical vari- 



ables in these solutions have power series expansions in terms of the harmonic functions 
AH X = H\ — 1, which all tend to zero as z-^oo lA . So in particular (noting the boundary 
conditions) we have the expansions 

= 1 + Atp + ... 
A x = k x (1-AH x + ...) 

h x = 1 + Ah x + ... (5.12) 

where the dots stand for higher than first order terms in AH x s. Using these in ( |5.3|) we 
find the expansion for V x (z) (to first order) as 

V\ = T x (^Ah\ — a x A(p + AH\) + ... (5.13) 

We identify the three leading order terms of this expansion as the gravitational, dilatonic, 
and (d\ + l)-form contributions to the potential energy of a distant (d x — l)-brane, which 
is parallel to the corresponding distribution within the system, that is 

V GX = l -T x Ah x , V DX = —T x a x A<p , V FX = T X AH X (5.14) 

Such an identification is possible, as each of these terms involves only the variations of 
the corresponding field variables and not mixings among themselves. As an example we 
consider the case of the orthogonal system di (1 d 2 = 5 and a (dx — l)-brane probe. Then 
from the solutions (|1.29| ) we find 



Ahi = -^AH 1 + 4(^-5)AH 2 

A(p = nlaxAHx + K\a 2 AH 2 (5.15) 

Specializing to the case of localized distributions along the transverse space at the points 
z x (A = 1, 2), we have 

AH x (z) = (5.16) 
\z - z x \ d 

where c x s are constants proportional to the densities of the corresponding longitudinal 
distributions. The potentials ( |5.14| ) for e.g. a (c?i — l)-brane probe at point z a thus become 

2 did x Tf 2l d x d 2 5 T{T 2 

V G i{z) = -a^—--- rf + C2K 2 ( 



D Iz-z^ zy 2D 2 ; \ z -z 2 \ d 

14 In general /S.H x s are written as linear combinations of h x s defined in (3.33). But for the orthogonal 
case considered here, they actually coincide. 
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xt / \ 2 a l T l 2 «l«2TiT 2 

V D i(z) = -CiK v 71-C 2 K 2 - -y 

\Z — Z\\ a \Z — Z2\ 

V Fl {z) = Cl ^L (5.17) 

\Z — Zi\ a 

where the first (second) term in each expression represents the interaction potential with 
the (di — l)-branes ((d 2 — l)-branes) distribution. As these formulas indicate, the grav- 
itational and dilatonic forces between two same-type (parallel) branes are always attrac- 
tive, while their form-field force is restrictly repulsive. In contrast, both the gravita- 
tional and dilatonic forces between two (orthogonal) different-type branes may be ei- 
ther attractive or repulsive depending on the values of the various dimensions and cou- 
plings, but there are no form-field forces between themselves, as is expected. When the 
(di — 1, di — l)-branes system is in a marginal (BPS) state, the internal branes as well as 
the brane probe do not feel any total force. So at an arbitrary point (z a ) we must have 
Vi(z) = Vqi(z) + Vdi(z) + Vfi(z) = 0, which according to ( |5.17| ) requires 

2 d\di 2 2,i n 

, d-\ do o . , . 

{ H~2 ) ~ aia2 = ° (5 ' 18) 

Interchanging the roles of {d\ — 1) and (d 2 — l)-branes in ( ^.18| ), we obtain the three 
consistency conditions ( 1.21[ ). This verifies explicitly our assertion in section 1 that the 



consistency conditions are nothings but the requirements of the no-force conditions be- 
tween different pairs of branes. 



6 Masses and charges 

The important physical parameters of a brane configuration (appearing e.g. in the black 
hole applications) are its mass and various charges. The total mass is the sum of the 
(positive) rest masses of the constituent branes, and the (negative) energy of the binding 
forces, both contained in the ADM mass of the brane system. For a marginal configura- 
tion, the total interaction energy vanishes and one expects that the ADM mass to be the 
sum of the constituents masses. The general ADM mass formula is derived by linearizing 
the Einstein equation in the (flat) background of the asymptotic metric | 4l| . For a system 



of distributed branes with a metric of the form (|3.1| ), the general formula (in the units 
with 167tGd = 1) reduces to 

M = = - f d i+2 zd 2 m(z) (6.1) 

Vd-l J 
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where Vd-\ is the volume of the distribution subspace, and ^(z) is 

= -f'hijiz) + h 00 (z) + (d+ l)e 2C(2) + const. (6.2) 

In this formula, 7^ = h t: > (z)\ z ^ 00 . Using the Stokes theorem, the expression for ADM 
mass becomes 

M = -f , d i+1 nr i+1 d r y(r,6 & ) = -^ +1 ( r " +1 <W™ (6.3) 

where (r, 9 a ) denotes a set of polar coordinates on the transverse space, and £^+1 stands 
for the unit (d + l)-sphere's area. From this formula a finite value for M. is obtained, 
only if \I/ at infinity behaves as \ (i.e. as an harmonic function). In fact by expressing 
the metric in terms of the harmonic functions h\ as in ( [3.33j ), and expanding in powers 
of h\s we obtain 

*(h x ) =J2^ + °( h D ( 6 - 4 ) 

A 

where fix's are constants depending on the specifications of the system under considera- 
tion. Using this in ( |6.3p , and reversing the route of Stokes theorem from ( |6.3j ) to ( |6.1D , 
yields 

M = -J^A*a / d S+2 zd 2 h x = J2 Hx I d i+2 z Px {z) (6.5) 

A J A J 

where we have used the Poisson equations d 2 hx = — px within the distribution region of 
the transverse space. Note that, despite using only the 0{hx) terms of (|6.4j) , this is an 
exact expression for Ai, since higher order terms of ^ do not contribute to the surface 
integral in (|6.3|). 



The charges are determined by types of the existing branes. In this paper we deal 
with cases where only electric charges present. The total electric charge, corresponding 
to a (d r + l)-form field strength jF r , is defined via a <i r -form conserved current J T as 

Q r = j *J r =f e 2a ^*r (6.6) 

where V dr+2 is any (d r + 2)-dimensional hyperplane intersecting all the existing (d r — 1)- 
branes at points and S dr+1 is a (d r + l)-sphere surrounding these points. The * in ( |6.6| ) 
denotes the Hodge dual in the curved background. By taking the radius of S dr+1 to 
infinity, the fields in (|6.6|) are replaced by their asymptotic values and it reduces to 

Q r = L*K,= I (6-7) 

where now * denotes the Hodge dual with respect to the flat background, and we have 
set as in the previous sections y?oo = 0. Using the ^4 r 's equation of motion at infinity (see 
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section 3) 

d*FL = *JL= E «aPaW*ca (6-8) 

d\=d r 

eq. ( |6.7| ) becomes 

Q r = = E *a / ^zpaW (6.9) 

where Vd-d r denotes the volume of the subspace transverse to [d r — l)-brane's world- 
volume and parallel to the distribution subspace. This shows that the charges of the 
same-type branes are additive even if they are not parallel, and n x p x (z) measures the 
charge density of the X th distribution per its unit transverse volume. We examine in the 
following the general formulas ( |6.9| ) and ( |6.5| ) for the orthogonal solution of section 2. 
Similar results hold for the non-orthogonal solution of section 4. 

The case of N orthogonal branes 

In this case from the ansatz (|2.2|) , and using ( |1.10|) and (1.12), we have 



= k x MH x A (A cfc*) (6.10) 



where the * denotes the Hodge dual in a (d+2)-dimensional Euclidean space (an irrelevant 
overall sign concerning dualization have been omitted). Using this in (|6.6|) , we obtain in 
agreement with Q6.9Q that 

Qa = «a / - *dH x = K X f d d+2 zp x (z) (6.11) 
Js^ 1 Jv d + 2 

where the Stokes theorem with and Poisson equation have been used. 
For the mass formula, using (|6.2|) with the solution (|2.12f) , we obtain 



* = E ^a 2 A# a + 0(AH 2 ) (6.12) 



A 



which shows that p x = k x 2 in this case. So by ( |6.5| ) we have 

M = Y J ^ 2 I d" +2 zp x (z) = E I^aQaI (6.13) 

A A 

That is the total mass M. equals the sum of the constituent masses M. x = \k x Q x \, which 
is an indication of marginality. 

Conclusion 

Starting from a reduced Lagrangian reformulation of the problem of orthogonal brane so- 
lutions, we arrived at a set of (linear) constraints, which was shown to consistently solve 
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the corresponding (non-linear) field equations. The requirement of consistency between 
these two sets of equations, led to a set of algebraic constraints containing all the physical 
information characterizing the marginal orthogonal solutions. These include the mass to 
charge ratios of the constituent branes and their suitable intersection rules. Although in 
the realistic supergravities this lead to the BPS saturated solutions with extremal (su- 
per) p-branes as the building blocks, extensions to the black solutions with non-extremal 



(black) p-branes are also possible using suitable deformation functions [22]. By intro- 
ducing a general formulation for handling arbitrary geometries of the intersecting branes 
with uniform 'longitudinal' distributions, a very general expression for the associated 
form-potentials in terms of the metric and dilaton field was derived. It was shown that 
the equations of motion of the reduced theory, can be translated to the 'forced-geodesic' 
equations describing a surface in the fields 'configuration space' . The conditions for the 
integrability of these equations are found to coincide with the constraints obtained earlier 
f34| . Essentially this type of formulation is not restricted to the case of the marginal 
solutions, as far as the number of the independent harmonic functions is not restricted to 
that of the density functions. The distributions densities may be so correlated to result 
in the dependent harmonic functions. Consequently the constraints of the marginal so- 
lutions are not valid for the non-marginal solutions. As a result the suitable intersection 
rules will be different from those of the marginal solutions. We hope that the formulation 
of this paper (with suitable changes) to be applicable for classifying these non-marginal 
solutions as well (see however P8fl). Finally we showed that how applying the ideas of the 
if-surface and null geodesic surface lead to the solutions for a system with two similar 
branes at SU(2) angles. 

Appendix: 

A. First order RL's for gravity 

The standard Einstein-Hilbert action is written in terms of the 2 nd order Lagrangian: 

C EH = V=gR(g,dg,d 2 g) (A.l) 

However, that the Einstein equation itself is of second order, shows that Ceh must be 
'equivalent' to a 1 st order Lagrangian. To see this explicitly, we apply the formula for 
Ricci tensor: 

1 

Z 9 

where the Christoffel symbols, T P MN 's are 1 st order quantities. Using this formula, the 
expression for y/—gR becomes 

V=gR = g MN d P {^-gT p MN ) - ^T P T Q PM - ^d M d M (\n J=g) 



R MN = ^=d P (^9T p MN ) - T P QM T Q PM - d M d N (ln y/=g) (A.2) 
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-gdpg 



MN r P 
1 MN 



-gT PQM T QPM + d M (ln^=g~)d N ( 
'N(V^99 MN d M ln^) 



-</</' ) 

(A.3) 



where the 2 nd order terms in the first line, have been appeared as total derivatives in the 
third line. After further simplification using the formulas 



9 



AfV-pP 

r 



°p9 l mn 

MN-rP 



QPR 



MN 



-2T PQR r 



equation (|A.3|) reduces to 



^R= C G (g } dg) - d M A M (g,dg) 



where we have defied 
C G 

A M : 



- T PMN T( 



(mn)p + d M (ln y/^g)d N ( 
^g-(d N g MN + 2g MN d N ln^9) 



~99 MN ) 



(A.4) 



(A.5) 



(A.6) 
(A.7) 



Formulas (|A.5|) to (|A.7|) can serve as simplifying explicit formulas for practical calculations 
of the Ricci scalar. However, since the surface term in (|A.5|) does not contribute to the 
equations of motion, we can identify the I s * order part, Cq, as the gravitational field's 
Lagrangian.[^] 



Application to semi-homogeneous spacetimes 

We have constructed examples of such spacetimes using the brane distributions through- 
out sections 1 to 4 of this paper. The generic form of the metric tensor for such spacetimes 
can be written as (see ( |3.1| ) and the related descriptions) 



9mn 



hij(z) 
5 ab e^ 



From this metric we have 



where h = \det(hi 



In = 1/2 In h + (d + 2)C 
The non-vanishing Christoffel symbols are 



aij 



l/2d a hi 



e 2C (S ab d c C - 5 bc d a C + 5 ca d b C) 



(A.8) 



(A.9) 



(A.10) 



15 It must be noted however that Cg despite Ceh is not transforming as a scalar density under coor- 
dinate transformations. 
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The expression for Cq in ( |A.7 ) thus becomes 



C G = h 1 l 2 e dC [l/Ad a h ij d a h i] + l/A(d a \nhf + (d + l)(d a \nh)d a C + d{d + l)(dC) 2 } (A.ll) 

After simple manipulations, using the matrix notation, this formula is rewritten as 

C G = e 2G [-l/ATrihr l db) 2 + l/4Tr 2 (h _1 ah) + (d+l)dC.Tr(br l dh) + d(d+l)(dC) 2 } 

(A.12) 

Here we have dropped the (contracted) transverse space indices, taken traces over 
and defined G as 

2G = l/2\nh + dC (A. 13) 

A simplifying trick for calculating Lq 

Finding a closed form for Lq m the specific problems (given the ansatz for Ay) by ( |A.12|) 
needs to a closed form of A y ", which in many situations can not be found easily. Fortunately 
a shortcut exists by means of which Lq can be computed without really inverting Ay. 
All that is needed, is to calculate the determinant: h = det(hij), and take its 1 st and 2 nd 
variations. This originates from the formula 

h " = ^ < A - 14 > 

dhij 

using which the two traces in ( |A.12| ), take simple expressions as 



- Tr(h-'dhf = ^^dh ir dh kl = (d 2 lnh) d2h ^ (A.15) 



Tr(h- L dh) = K '- dhij = d(lnh) 

Oflij 



That is the two traces are calculated by expressing In h as a function of (Ay), taking its 
I s * and 2 nd order variations, and neglecting the 2 nd order variation of 0. Finally these 
two expressions are combined into the equation 

Tr 2 (h-^h) - rr(h _1 9h) 2 = ( ^ J (A.16) 
This is proved to be a very useful formula for the sake of practical calculations. 

Application to the system d\ fl d 2 = 5 at angles 

For this system Ay- is defined by ( |4.1|) . Computing the determinant of this matrix we 
obtain 

h = /( g )e 2 ( 5B °+5iBi+&B 2) (A. 17) 



16 



When calculating the second variation, we are allowed to replace hij's with another set of variables 
in terms of which hij's are linear. 
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where 

f(q) = \det(l-q 2 11 T )\ (A.18) 
Defining the four independent variables q A linear in hi/s as 

q A = (e 2B \e B \e 2B \qe B ^) (A. 19) 

and using the above tricks, with d 2 q A = in ( A.16|) , we obtain 



l/2Tr(h- 1 9h) = 5dB + 8 l dB 1 + 5 2 dB 2 + ^-dq 

2 J 

l/A[Tr 2 ihr l dh) -Trih^dh) 2 } = 

5(5 - l)(dB ) 2 + M5, - 1) + ^(dB,) 2 + [5 2 (5 2 - 1) + ^j](dB 2 ) 2 + 

af 

255 1 dB .dB 1 + 255 2 dB .dB 2 + 2(5 1 5 2 - jj)dB v dB 2 + 

f jdq.[5dB + (5 1 - l/2)dB 1 + (5 2 - l/2)dB 2 ] + f -^{dqf (A.20) 



Using these in (|A.12| ) and including the dilaton term, a final expression for Lq as in (f4.4|) 
results. 

B. Analysis of the Diophantine equation for intersections 

As we have seen in section 1 , possible marginal intersections of super p-branes are governed 
by a Diophantine equation written as ( |1.28| ) or 

(2D - d 1 d 1 )(2D - d 2 d 2 ) = (d x d 2 - 5D) 2 (B.l) 

There are at least two means for classifying the solutions of this equation. Given the 
spacetime dimension D, we can specify: 

1) the number of common directions (5 — 1) or 2) the number of angles m 
and look for the possible dimensions of the branes (d± — l,d 2 — 1). We first present 
the method of analysis for arbitrary D's below, and at the end summarize the results for 
interesting dimensions D — 4, 6, 10, 11. But before, two simple cases may be distinguished: 
a) The same-type branes 

This corresponds to: d\ = d 2 & «i = a 2 so that a\a 2 > and ( |1.21| ) implies: 2D — d\d\ = 
+ (dj - 5D) , i.e. 

5 = d l -2 (B.2) 

This means that for two branes of the 'same type' to marginally bind, all their angles 
except two of them must be vanishing (i.e. m = 2). A result which is in section 4 of this 
paper. 
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b) The self-dual pair of branes 

This case corresponds to: d\ = d 2 & oc\ = —a 2 so that a\a 2 < and (|1.21|) implies: 
2D — d\d\ = —{d\di — 5D) , i.e. 

5 = 2 (B.3) 

This means that for a 'self-dual' pair of branes to marginally bind, they must be inter- 
secting over a string. 

General restrictions 

Not all the solutions of (P-l| ) are physically acceptable, as we have two sets of restrictions: 
First, the definitions of dimensions and the condition of asymptotic flatness require 

0<5<d x <D ; A = l,2 

0<d 1 + d 2 -5 <D (B.4) 
Second, the reality of («i, a 2 ) and the restriction on their relative sign require 

d\d\ <2D ; A = 1,2 

sgn(did 2 — SD) = sgn(aia 2 ) (B-5) 

The method of classification by 5 

Noting the symmetry of ( |B.1| ) relative to (di,d 2 ), we introduce the new variables 

r = did 2 , s = d 1 + g? 2 (B-6) 
in terms of which ( |B.1|) is written as 

r = 2s + P + — %- (B.7) 
s-Q 

where (P,Q,R) are integers defined as 

P = 4(5-2) , Q = D + 2(5 -2) , R = (8 - D)(5 - 2) 2 (B.8) 

Therefore solving (p.l| ) for (di,d 2 ) (assuming (5, D) are given) is equivalent to solving 
( |B.7| ) for (r, s) generally, and selecting then those solutions for which (d\, d 2 ) are integers, 
i.e. 

\ (2q + l) 2 , SG 2Z + 1 W 7 V ; 

On the other hand ( |B.7| ) by itself has a finite set of solutions for (r, s) , for given values of 
(P, Q, R), which may be obtained by demanding that the integer (s — Q) enumerates the 
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integer R ^ 0. 



The R = cases: 

These include: D — 10 and/or 5 = 2. In both cases each of the relations: s = Q 



r = 2s + P solves the equation (|B.7|) . The corresponding solutions for (di,d 2 ) then are 
given by 



D = 10 : \ or 



d 1 + d 2 = 2(5 + 2) 
(d 1 -2){d 2 -2) = A{5-\) 



(B.10) 



5 = 2 : 



or 



d\ + d 2 = D 



(di,d 2 ) = (4,4) J (3,6) 



(B.H) 



The method of classification by m 

For simplicity we may assume: g?i < d 2 and thus 5 = d\ — m. So noting (|B.6|) and (|B.8|) 
we can write ( |B.7[ ) as 

(8-L>)(di - m-2) 2 



did 2 = 2(di + d 2 ) + 4(di - m - 2) + , 

(a 2 — cti) — + 2(m + 2J 

Defining the 1 to 1 map: (di, d 2 ) i— > (x, y) in the domain of Z by 

a; = di - (m + 2) , y = (d 2 — d{) — D + 2(m + 2) 
equation ( B.12 ) in (x, ?/) variables transforms to 



x 2 + yx + (my — k — 



I 







y-n 



where (n, fc, /) are integers defined by 

- D , k = (m — 2) 2 , Z = rafc 



n 



;b.i2) 



(B.13) 



(B.14) 



(B.15) 



Solving ([B.1D for (d 1; d 2 ) (assuming (D, di — S) as given) is equivalent to solving ( p3.14j) for 
(x, y) which is much easier. In fact equation ( |B.14j ) (as ( |B.7| )) has a finite set of solutions 
for given values of (m,n,k,l), which can be found easily by demanding that the integer 
(y — n) must enumerate the integer I ^ (for I = see below), and further that 



y 2 — A[my — k 



I 



y-n 



(2g) 2 
(2g + l) 2 



s G 2Z 
s 6 2Z + 1 



(qeZ) 



(B.16) 
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so that ( |B.14| ) has integer solutions for x. 



The I = cases: 

By ( p,15|) these include: D = 10 and/or m = 2. Again (like the R = case) the solution 

t(m-l) 



in both cases lie on two branches, i.e.: y = n , y = —x + m — 4 ^ 1 K The corresponding 



solutions for d 2 ) then are given by 

!d 2 - dx = 4 - 2m 
or (B.17) 
(di-2)(d 2 -6) = 4 -4m 

or (B.18) 
(di, d 2 ) = (3, D - 8), (4, D - 6), (6, £> - 5) 

where the condition c?i < d 2 in the second class of (|B.18| ) requires that D > 11,10,11 
respectively. Note that the two classes: (B.17) & ( |B.1U ) are the same. The condition 



d\ < d 2 restricts the first class of the ( |B.17| ) solutions to those with m = 0,1,2. In 



particular for two parallel branes (m = 0) or two non-parallel branes with only one angle 
(m = 1) in this class, we obtain: d 2 — d\ =4 or 2 respectively, in agreement with the 
results of [[57] for D-branes. In fact noting that the dilaton coupling for a D- {d — l)-brane 
is a(d) = ±(4 — d)/4, with +(— ) sign for branes with electric (magnetic) RR charges, we 
conclude that the first (second) class of solutions in ( |B.17[ ) in the case of two D-branes, 
describes bound states of two D-branes of the same (opposite) 'electromagnetic' type. 



Special cases of D = 4, 6, 10, 11 

We present here the summary of the above classifications of solutions for D = 4, 6, 10, 11. 
In this summary we relax the restrictions: 5 > , d < D to the extent that: 5 > , 
d < D; so that the solutions include 'instanton-like' objects as well as the 'non- asymptotic 
flat' configurations. The last column in each table idicates that the two branes are of the 
same or different electric/magnetic type according to ai<y 2 to be positve or negative. 



5 


(di,d 2 ) 


m 


sgn(aia 2 ) 





(2,2) 


2 


+ 


table (2) D = 4 


S 


(di,d 2 ) 


m 


sgn(aia 2 ) 





(2,2) 


2 


+ 


1 


(3,3) 


2 


+ 


2 


(2,2) 
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table (3) D = 6 



5 


(di,d 2 ) 


m 


sgn(aia2) 





(0,4), (1,3), (2,2), (1,6) 


0,1,2,1 


0,+,+,+ 


1 


(1,5), (3,3), (2,2), (2,3),..., (2,9) 


0,2,1,1,...,1 


-,+,-,-,0,+,...,+ 


2 


(2, 6), (3, 5), (4, 4), (3, 6) 


0,1,2,1 


- - o,+ 


3 


(3, 7), (4, 6), (5, 5) 


0,1,2 


-,o,+ 


4 


(4, 8), (5, 7), (6, 6), (5, 6) 


0,1,2,1 


0,+,+,- 


5 


(5, 9), (6, 8), (7, 7), (6, 6) 


0,1,2,1 


+,+,+,- 


6 


(6,7) 








table (4) D = 10 



5 


(d 1 ,d 2 ) 


m 


sgn(a:ia 2 ) 





(0,3), (2, 2), (0,6) 


0,2,0 


0,+,0 


1 


(3,3) 


2 





2 


(2, 7), (3, 6) 


0,1 


- o 


3 


(3,9) 








4 


(6,6) 


2 





5 


(7,7) 


2 


+ 


6 


(6, 9), (8, 8) 


0,2 


o,+ 



table (5) D = 11 



C. Definition of the angles between two branes 

We have seen at the beginning of section 4 that the 'relative' orientation of a pair of 
branes can be described in terms of the parameters 7 mm / entering the metric as in ( |4.1|) . 
The asymptotic form of this metric, restricted to the subspace of coordinates {vx,y%') 
then is written as q] 

da 2 = (dy?) 2 + (dyf) 2 + 2 lmm ,dy™dy? (C.l) 
Therefore 7 mm ''s are related to the angles between coordinates as 

7mm' = COs(d m , 9 m ')L-»oo ( C - 2 ) 

Clearly the (^i^) numbers 7 mm ' depend on the choice of the (Cartesian) coordinate sys- 
tem: (y™, y™'), transforming as components of a rank (1, 1) tensor under the rotations: 
SO{5\) x 50(^2) of these coordinates. Therefore we can not identify (j mm i) as the set 
of independent parameters needed to describe the relative orientation of the two branes, 

17 Refer to table (1), section 1, for the definitions of coordinates and subspaces. 
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since they are related through these rotations. To do this, we need to specify the maximal 
set of SO(5i) x SO {62)— invariant quantities. We call such invariant parameters as the 
'geometric' or 'intrinsic' angles of the two branes. To give a simple description of these 
angles, we use our intuitions in 3-dimensional Euclidean geometry^. We first take a pair 
of arbitrary unit vectors (ni,n 2 ) within the (61,62) subspaces respectively as 

m = u?d m ; «) 2 = i 

n 2 = wfflU' ; «') 2 = i (C3) 

The angle between (ni, n 2 ) according to the metric (|C.1| ) then is defined by 

cosO = ni.n 2 = ^ mm >uj'1 l uj 2 n ' (C.4) 

We now define the 'geometric' angles between the two branes to be the 'non-trivial' 
extremums of the quantity A = cosO. By 'non-trivial' here, we mean those extremums 
which are not zero identically, and not related together by a change of signs. These 
extremums are obtained by extremizing the function 

at ,\ ~ .Ti ,m' ^1 / ,ni\2 "^2 / m' \2 I /~i c\ 

b(u) = Jmm'Wl U 2 - —(U l ) - y(u; 2 ) (C.5) 

subject to the constraints in ( |C3| ) with (Ai,A 2 ) as the Lagrange multipliers. For the 
extremum points we have 

Ai = A 2 = A = cosO (C.6) 
where A is obtained from the eigen-value secular equation 

-(1 + A)1]=0 (C.7) 



where 7 = ( T ) is a (61 + 6 2 ) x (6\ + 6 2 ) matrix representing the metric tensor 



1 

1 , 

of ( |C.1| ). So (|C.7| ) is written as a polynomial equation in A of degree (6\ + 6 2 ), whose 
'non-trivial' roots {A,.},.^...^ give, via A r = cos8 r , the geometric angles {8 r }. Not that 
by Hermiticity of 7 all the A r 's are real. Also by positive-definite ness of da 2 we can show 
that always — 1 < A r < 1. Further the secular equation (|C.7|) can be written as 



\ S2 ~ 5l det(jj T - \ 2 ) = (C.8) 

which (assuming 61 < 62) implies that the roots of (|C.7| ) consist of a set of (62 — 6\) zeros 
and 61 pairs of opposite numbers (this had been expected, since reversing the direction of 
ni or n 2 in ( |C.4|) changes the sign of cosO but preserves its extremum property). So the 

18 You may consider e.g. a line and a plane at angle in the 3-space. 
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number of # r 's equals the degree of the determinant in ( p.8|) as a polynomial function of 
A 2 , which is in general 

m = Min{8 1 ,8 2 } (C.9) 

This number is in fact the (maximum) number of the successive rotations, required for 
bringing the smaller in dimension brane from a parallel to an angled status relative to the 
other brane. 



D. Derivation of useful i7-surface identities 

The model for the form-field sector of the RL ( |3.40| ) can be simulated by a simplified 
model in the discrete mechanics as 



(D.l) 



where (x a ) and (y l ) are two sets of dynamical variables and gij{x) is any invertible 'metric 
tensor' . The equations of motion for y l, s are obviously integrable and yield 



(D.2) 



where c/s are some integration constants. Eliminating y l, s in (|D.1|) by this equation, we 
obtain 

L[x] = ~c iCj g^{x) (D.3) 



Applying equation (]D.2| ) again, it is easy to prove the 'on-shell identity' : 



_d_ 

dx° 



L[x,y] 



d 



dx c 



L\x] 



(D.4) 



Changing the role of (x a , y\ 7^, L) to (^f> a , Afy, e 2(Xrtp (h l i), Cf^J in this equation, and using 
the definitions (|3.44j ) and (|3.50|) , we obtain the 'H-surface identity' : 



Ux\', a ((p,H) = -U X \> 



(D.5) 



Differentiating the if-surface identity (|3.49|) relative to H\u and using ( p.5|) , we obtain 
another identity: 

U aa',a" (0) H) = 2d xll u xy ((j ) ) = 2d xll U xy ((f), H) (D.6) 
Using this result in the equation Q3.47 ), after simple manipulations we obtain 



dx>U: 



X'^XX" 



dx"Uxx' 



(D.7) 



Noting the symmetry of U\\> in its two indices, this equation implies the existence of a 
function U(H), such that 

U xy (<f>(H),H) = d x d xl U(H) (D.8) 
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This result could be seen in a (somehow) different manner. Let's define 




(0)d x AUH)d x ,AUH) 



(D.9) 



so that U- 



E"=i U r xy by (H§. From (ggg) it is evident that 



U r xx ,(4>(H),H) = d x (c^A r (i) ) = 



(D.10) 



This implies the existence a function U r (H), in terms of which we can write 



4%) = dxU r 
U r xx , = d x d x ,U r 



(D.ll) 
(D.12) 



which proves (|D.8| ) again. 
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